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This  is  a  collection  of  lectures  delivered  at  the  Universtdade 
Federal  de  Pernambuco,  Recife,  Brazil,  between  2  November  1981  and  2 
December  1981*  The  visit  to  Recife  was  a  response  to  an  invitation  of 
Professor  Gilberto  F.  de  S&  of  the  Physics  Department.  In  the  prepara¬ 
tion  of  these  notes  I  made  many  requests  of  the  research  workers  still 
at  HDL  for  earlier  results  and  recollections  of  our  early  work.  Among 
those  most  frequently  consulted  were  Donald  Wortman  and  Nick  Karayianis 
and  I  wish  to  thank  them  for  their  prompt  responses  and  their  valuable 
suggestions.  The  later  work  and  a  critical  reading  of  the  entire  rough 
notes  of  the  lectures  were  done  by  my  coworker  Richard  P.  Leavitt.  To 
him  I  owe  a  special  thanks  for  numerous  helpful  suggestions.  A  number 
of  suggestions  from  my  Brazilian  colleagues  helped  make  the  lectures 
more  clear.  Particular  among  these  were  Professor  Oscar  Malta  and 
Professor  Alfredo  A.  da  Gama  both  of  whom  I  wish  to  thank  for  their 
help.  1  would  also  like  to  thank  Evandro  J.  T.  de  Araujo  Gouveia  for 
his  suggestions  and  finding  errors  in  the  notes.  In  al.l  of  the  efforts 
to  make  the  lectures  more  complete  I  owe  special  thanks  to  Professor 
Gilberto  F.  de  SA.  Lastly  I  would  like  to  thank  the  entire  physics 
department  at  the  University  of  Recife  for  making  our  stay  in  Recife  a 
warm  and  pleasant  experience.  Both  my  wife  and  myself  have  never  before 
experienced  such  a  pleasant  display  of  hospitality.  To  all  the  members 
of  the  physics  department  my  wife  and  I  say  "muinto  obrigado." 
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INTRODUCTION 


The  work  at  the  U.S.  Army  Harry  Diamond  Laboratories  (HDL)  on  the 
experiment  and  theory  of  rare-earth  ions  in  transparent  host  materials 
was  begun  in  the  early  1960's.  This  was  after  workers  at  Bell  Telephone 
Laboratories  (Johnson  et  al,  1962)*  demonstrated  the  first  continuous 
operating  laser  using  neodymium  (Nd)  in  calcium  tungstate  (CaW04).  It 
was  decided  at  that  time  that  the  HDL  group  would  concentrate  on 
NdsCaWO^.  The  group  was  limited  in  number  at  any  given  time,  but  over 
the  span  of  years  many  research  workers  were  participants.  A  list  of 
publications  pertinent  to  the  theoretical  analysis  and  experimental  work 
on  triply  ionized  rare-earth  ions  in  solids  is  given  in  the  annotated 
bibliography. 

The  plan  that  the  group  was  to  follow  during  the  first  two  years  and 
still  follows  is  the  following:  (1)  As  many  as  possible  of  the  rare- 
earth  ions  were  to  be  doped  into  the  host  material  CaW04,  and  the  opti¬ 
cal  data  were  to  be  carefully  recorded.  (2)  These  data  were  to  be 
analyzed  by  any  techniques  existing  in  the  literature  or  developed  at 
HDL  to  obtain  phenomenological  crystal-field  parameters.  (3)  A  theo¬ 
retical  derivation  of  the  crystal-field  parameters  was  to  be  developed 
which  gave  not  only  the  parameters  which  were  obtained  from  the  experi¬ 
mental  data,  but  also  a  set  of  odd-k  crystal-field  components.  (4)  The 
odd-k  crystal-field  components  were  to  be  used  to  calculate  the  electric 
dipole  transition  probabilities  using  the  theory  of  Judd  (1962)  and 
Ofelt  (1962). 

As  of  1970,  all  these  objectives  had  been  met  for  CaW04  (and  conse¬ 
quently  to  a  lesser  extent  for  any  solid  where  the  rare-earth  ion  occu¬ 
pies  a  site  of  S4  symmetry).  During  this  time,  a  number  of  new  theories 
as  well  as  significant  modifications  of  existing  theories  were  devel¬ 
oped.  The  most  important  theory  thus  developed  was  the  HDL  three- 
parameter  theory  of  crystal  fields.  This  theory  was  a  blending  of  the 
much  older  point-charge  model,  including  the  effects  of  covalency  and 
wave-function  expansion  and  screening  of  the  field  of  the  point  charge. 
These  effects  were  encompassed  in  the  three  parameters,  and  a  consistent 
set  of  these  parameters  was  found  for  the  triply  ionized  rare-earth 
series  in  the  host  CaW04.  The  phenomenological  crystal-field  parameters 
used  in  the  development  of  the  theory  were  those  obtained  by  using  an 
effective  spin-orbit  Hamiltonian  developed  by  Nick  Karayianis  at  HDL. 
The  optical  data  were  taken  by  Donald  E.  Wortman  and  Ruben  T.  flarrar  of 
HDL. 


*Because  of  the  large  number  of  literature  citations ,  these  are  given 
in  abbreviated  form  only,  not  footnoted  on  each  page.  Complete  liter¬ 
ature  references  are  given  alphabetically  by  author  in  the  Literature 
Cited  section. 


9 


Since  1970,  a  number  of  new  additions  have  been  made  in  the  theory, 
and  the  ability  to  analyze  the  optical  data  has  been  achieved.  These 
include : 

(1)  The  development  of  lattice  sums  (point-charge  model)  for  any 
crystalline  solid  (230  different  space  groups). 

(2)  The  self-consistent  inclusion  of  the  dipole  fields  for  any 
crystalline  solid. 

(3)  The  calculation  of  the  phenomenological  crystal-field  param¬ 
eters  from  experimental  data  for  any  solid  having  point  sym¬ 
metry  higher  than  C1  or  . 

(4)  The  calculation  of  transition  probabilities,  Zeeman  splitting 
factors,  Judd-Ofelt  intensity  parameters,  branching  ratios,  and 
lifetimes  for  any  point  symmetry  higher  than  or  . 

Of  the  recent  work  at  HDL  done  on  rare  earths,  the  most  significant 
is  the  work  done  on  rare-earth  ions  in  lanthanum  trifluoride  (Ln:LaF.j, 
Ln=rare-earth  ion)  (C2  symmetry)  and  reported  in  two  papers.  Marrison 
and  Leavitt  (1979)  reported  an  analysis  of  the  low-temperature  optical 
data  which  had  been  previously  reported  by  workers  in  other  laboratories 
(very  much  by  W.  T.  Carnall  of  Argonne  National  Laboratory).  In  our 
analysis,  the  point-charge  lattice  sums  were  used  to  obtain  starting 
crystal-field  parameters  in  a  least-squares  fitting.  Consequently,  sets 
of  crystal-field  parameters  for  almost  all  the  rare-earth  ions  were 
obtained  which  were  consistent  through  the  entire  series.  In  the  second 
paper  (Leavitt  and  Murrison,  1980),  the  results  of  the  lattice  sum  for 
the  odd-k  crystal-field  components  were  used  to  calculate  the  Judd-Ofelt 
intensity  parameters  and  lifetimes  for  the  entire  rare-earth  series. 
Wherever  possible,  these  were  compared  to  experimental  data,  and  the 
agreement  was  found  to  be  very  good. 

At  present,  with  J.  B.  Gruber  and  N.  C.  Chang,  we  are  analyzing  data 
on  triply  ionized  rare-earth  ions  in  yttrium  oxide  (Y^^).  This  mate¬ 
rial  has  two  types  of  sites  that  the  rare-earth  ion  can  occupy,  C2  and 
C^.  A  paper  on  "The  Analysis  of  the  Spectra  of  Kramers  Ions  in  the  C2 
Sites"  has  been  published  (Chang,  Gruber,  Leavitt,  and  Morrison,  1982). 
A  second  paper,  "Optical  Spectra,  Energy  Levels,  and  Crystal-Field 
Analysis  of  Tripositive  Rare-Earth  Ions  in  Y20^.  Non-Kramers  Ions 

in  Cj  Sites"  has  been  published  (Leavitt,  Gruber,  Chang,  and  Morrison, 
1982).  We  are  now  writing  a  third  paper  on  the  rare-earth  ions  in  Y-O-j 
for  the  rare-earth  ions  in  the  sites.  This  latter  paper  also  will 

summarize  all  the  work  on  C2  sites. 

If  we  find  time  and  money,  we  intend  to  extend  our  computation 
ability  to  incorporate  more  of  the  free-ion  interactions.  Also,  we 
would  like .  to  incorporate  in  our  programs  many  of  the  theories  that  we 
have  developed  at  HDL  over  the  last  few  years.  The  inclusion  of  the 
self-consistent  dipole  contributions  to  our  theory  of  crystal  fields 
should  improve  both  the  energy-level  calculations  and  the  transition 
probabilities. 
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2.  ANGULAR  MOMENTUM  ALGEBRA 


2 . 1  Angular  Momentum  Operators 

In  classical  mechanics,  the  angular  momentum  of  a  particle  is 
defined  by 


I  = 


-► 

r 


x  p  • 


( D 


Actually,  we  should  specify  that  the  angular  momentum  so  defined  is 
about  a  particular  origin,  and  r  is  the  vector  distance  from  the  refer¬ 
ence  point  to  the  particle  with  momentum  p. 


If  we  use  the  commutation  relations 


[  x . ,  p  .  ]  =  rfiiS .  .  (  2 ) 

with  x^  =  x,  y,  or  z,  then  we  can  obtain  the  commutation  rules  for 
angular  momentum. 


\l  ,1  1  =  ifU  ,  [A  ,t  ]  =  itU  ,  and  ,1  1  =  idU  ,  (3) 
i  x  y J  z  1  y  z J  x  LzxJ  y 

which  are  the  basic  commutation  rules  for  the  Cartesian  components  of 
the  angular  momentum.  For  convenience  here  we  shall  drop  the  -fi  in  the 
commutation  relations.  This  does  not  mean  that  we  drop  -fi  throughout;  v?e 
restore  -ft  simply  by  writing  the  interactions  involving  the  angular 
momentum  so  that  the  -R  is  contained  in  the  constants.  As  an  example  of 
this,  we  consider  the  spin-orbit  Hamiltonian 


H 


2 


1 

2mzc2 


130  t  + 

-  T—  i’B 

r  3r 


(4) 


with  1  and  s  (the  spin  angular  momentum,  s,  we  will  discuss  later) 
containing  fi  as  in  equation  (3).  When  written  in  terms  of  unitless  t 
and  s,  we  have 


H 


2 


*z 

2m2c2 


1  3U  t  + 
-  c-  Jt*s 
r  9r 


(5) 


where  1  and  £  obey  the  commutation  rules  in  equation  ( 3 )  but  'R  =  1 . 

For  our  purposes  here,  we  want  to  use  the  spherical  represen¬ 
tation  of  I,  which  is  given  by 
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and  the  commutation  rules  are 


f  4+1 ' ] 


*+1  ' 


*-1  ' 


The  eigenfunctions  of  the  angular  momentum  are  the  spherical 


harmonics,  Y.  (9,<{>)  and 
Jem 


iz\ Jim>  »  m| £m>  , 


_  [  (  l+m)  (  l±m+1 )  ]  2  |£;m±1>  , 

1 1  /o 


where 


lm>  =  Y  (6f«|i) 


Frequently,  we  shall  use  the  unit  vector  r  to  indicate  the  argument 

of  Y.  ,  thus: 

Jtm 

When  the  Y„  are  wave  functions  such  as  in  equation  (8),  we  have 

Y£mU)  =  Um>  . 

The  wave  functions  have  the  property  that 
<£,m'l*m>  =  6U’6mm'  * 

Further,  we  shall  assume  that  the  spin  angular  momentum,  s,  obeys  the 
same  commutation  relations  as  equation  (7);  the  two-component  spinor 
wave  functions  are  represented  by  the  wave  function  |smg>,  so  that  the 
single  electron  wave  function  for  orbital  and  spin  angular  momentum  is 

|fcm£>|sms>  .  (9) 

The  wave  functions  given  by  equation  (9)  then  obey  the  following: 

l  |lm.>|sm  >  =  m„l£m.>jsm  >  , 

o  1  s  i  l  s  (10) 

( I) 2 1  lm.>  |  am  >  =  l( l+l ) | lm  > |sm  >  , 
xs  K  s 
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s  i  Jim. >  Ism  >  =*  m  |im.>|sm  > 
o  i  s  s  t,  s 


(s)2|f,m  >|sm  >  =  s ( s+1 )  |  Im  >  | sm  >  , 

X.  S  x>  s 

where,  of  course,  s  =  1/2.  A  further  property  of  the  spherical  har¬ 
monics  is 


A  A 

where  the  inversion  operator  is  Ir  =  -r,  a  property  that  will  be  used 
frequently  in  our  analysis.  While  many  of  the  interaction  terms  of  the 
Hamiltonian  were  derived  by  using  spherical  harmonics,  it  is  convenient 
to  introduce  the  ter.sor  operators 

cte|'>  *  /3S  w;>  <,2> 

Since  Y^lr)  =  (-DmYJlf_ra(r)  ,  we  have 

C**m<*)  -  *  <13> 

The  use  of  rather  than  Y^  in  the  interaction  terms  eliminates 
almost  all  the  factors  of  4ii,  An  example  of  this  is  the  coupling  rule 
for  spherical  harmonics  (Rose,  1957:  61): 


Y  v  =  Y  r(2k+1)(2n-*-1)-|I/2 

kq  nm  £  L  4ir(2*+i)  J  {14) 

x  <k(0)n(0)  |i(0)xk(q)n(m)  |jl(q+m)>YJl  q+m  , 

but 

C.C  *  I  <k(0)n(0)U(0)Xk(q)n(m)U(q+m)>C.  .  (15) 

kq  nm  r  i,q+m 


In  equations  (14)  and  (15)  all  the  tensor  operators  have  the  same  argu¬ 
ment  . 


2.2  Clebsch-Gordan  Coefficients 


For  our  purpose,  it  is  convenient  to  define  the  Clebsch-Gordan 
(C-G)  coefficients  as  the  coefficients  in  the  transformation  from  two 
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angular  momentum  spaces,  say,  1  and  s,  to  form  the  composite  space 
That  is, 

I  jm>  -  l  <£(  p)s(m-p)  |  j(m)  >  |Jtp>  |s,m~p>  ,  (16) 

U 

where  the  quantity  <  £(  p)s(m-p)  |  j  (m)  >  is  a  C-G  coefficient.  Since  we 
wish  to  have  ortnonormal  basis,  we  have 


since 


< j'm’ | jm> 


5.  .,6  , 

j  j '  mm' 

l  <fc (p )s(m-p ) | j (m) ><£ (p ' )s(m'-p ' ) | j ' (m' )> 

UU  ' 


(  17) 


< Jlp '  I  £ »j>  =  6  and  <s ,m' -p '  |  s ,m-p>  =  6  ,  ,  (18) 

.  Up'  m-p,m'-p' 

Thus,  we  have 

8  .  . ,  =  l  <£( p)s(m-p) | j(m)><£(  p)s(m-p)  |  j' (m)>  ,  (19) 

33  P 

an  important  and  very  useful  result.  If  we  assume  (correctly)  that  the 
same  coefficients  connect  the  j  space  to  the  1  and  s  spaces,  we  can 
obtain  another  condition  on  the  C-G  coefficients,  namely, 

6£r6ss'  =  E  <£(m£)s(ma)|  j(mJl+ms)><£’(mJl)s,(ms)|  j(m£+mg)>  .  (20) 

j 

Some  other  relations  among  C-G  coefficients  are 

<a(a)b(  p)  |c(  y)>  =  0  if  |ot|  >  a,  or  |Bl  >  b,  or  |y|  >  c 

and  if  y  *  a  +  8  • 

The  C-G  coefficients  vanish  unless  the  three  angular  momenta  obey  the 
triangle  condition,  or  |a-b|  <;  c  £  a  +  b  and  any  permutation  of  a,  b, 
or  c.  Further  properties  of  the  C-G  coefficients  are  given  in  appendix 
A. 


The  commutation  relations  for  the  spherical  components  of  the 
angular  momentum  of  a  single  electron  given  in  equations  (7)  and  (8)  can 
be  written  compactly  in  terms  of  C-G  coefficients  as 

[*p,*v]  -  /2  < 1 ( v) 1 ( w) I 1(p+v)>£y+v  (7a) 
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and 


[a  ,s  ]  -  /2  <1(v)1(u)!l(u+v)>s  .  (8a) 

1  P  VJ  y+V 

The  total  angular  orbital  momentum  operator  for  a  system  of  N  electrons 
is 

N 

L  -  l  l(i)  ,  (21) 

i 

and  the  total  spin  angular  momentum  operator  is 

N 

I  -  l  s(i)  .  (22) 

i 

The  spherical  components  of  these  operators  obey  the  same  commutation 
relations  as  equations  (7a)  and  (8a),  or 

[h  ,L  ]  =  /2  < 1 <  v) 1 ( y ) | 1(u+v)>L  (23) 

L  y  vJ  y+v 

and 

[s^sj  -  /2  <1(V)1(|i)|l(y+v)>S  . 

Also,  it  should  be  noted  that 

[lyiv(i>]  =  /2  <1(v)1(y)  |  1(y+v)>Zjj+v(i) 

and 

[Sy,sv(i)  ]  =  /2  <1(v)1(|i)|1(y+v)>s  <i)  .  (26) 

Consequently,  from  using  equations  (23)  through  (26),  we  have 

[lyCkq(i)]  =  /kTir+TT  <k(q)1(y)|k(q+u)>Ck^q+ij(i)  ,  (27) 

tw1’1  ‘  (wi,] 

with 

1  =  £  +  £  . 

i 

t 
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(28) 


(29) 


(24) 

(25) 


2*3  Wigner-Echart  Theorem 


The  Wigner-Echart  theorem  states  that  if  we  have  a  spherical 
tensor  Tkq  in  the  space  spanned  by  the  wave  functions  |JM>,  then  the 
matrix  elements  are 

<J’M*  |Tkg|JM>  =  <J(M)k(q)!j,{M,)XJ’ «Tk»J>  .  (30) 

The  projection  (q)  dependence  is  contained  in  the  C-G  coefficients,  and 

the  factors  <J'IT,  (J>  are  called  the  reduced  matrix  elements, 
k 

If  we  have  a  mixed  spherical  tensor,  rank  <,  projection  X,  in 
spin  space  and  rank  k,  projection  q,  in  orbital  space,  the  Wigner-Echart 
theorem  then  is 

<L,"Ls'Ms<|LVV 

-  <L(ML)k(q)|L'(M^)><s(Ms)lcU)|S'(M^)>  (31) 

x  <L'S*  IT*CklLS>  . 

Since  the  C-G  coefficient  is  purely  a  geometrical  factor,  all  the  phys¬ 
ics  is  contained  in  the  reduced  matrix  element.  The  Wigner-Echart 
theorem  allows  the  extraction  of  the  geometrical  factors  from  many 
complicated  matrix  elements;  it  also  serves  as  perhaps  the  main  motiva¬ 
tion  for  the  development  of  Racah  algebra  in  dealing  with  angular  momen¬ 
tum  states. 

Because  of  the  power  of  the  Wigner-Echart  theorem,  it  occurred 
to  Racah  to  cast  the  various  operators  representing  the  interactions  in 
terms  of  universal  quantities  that  could  be  tabulated  for  a  frequently 
used  many-particle  system.  Toward  this  end,  Racah  introduced  the  unit 
spherical  tensors  for  the  electronic  configuration  niN  which  we  define 
as 


<A’m'|ukq|fm>  -  <£(y)k(q)U(m')>6u, 


for  the  orbital  space  and 


<£ 'm' s 'mi |  v*k  |£msms> 

=  <£  (m)k(q)  jit  (m'  )><s(mgjK(X  )  |s(m^)  > 
x  6££*6ss’  ' 

for  orbital  and  spin  space. 


(32) 
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The  generalization  to  an  N-electron  system  is  simply 


'Jkq  ‘  j  V1’ 


A  simple  and  often  used  example  of  these  tensors  in  orbital  space  is  (we 
shall  omit  the  upper  limit  on  the  i  sum  in  the  remainder  of  the  discus¬ 
sion) 

l  V11  - 1 

»  (34) 

«  <  i  1C.  I  fc>U. 

k  kq 


where 


<AIICkH£>  -  <£(0)k(0)  U(0)> 


h  =  T  l  (i) 
y  |  y 


-  I  <tl£ll>u1  (i) 
i  14 

=  <HSLt  1>U,  ,  and 

ly 


<mu>  -  /*u+n  . 


An  example  of  a  tensor  in  a  mixed  spin  and  orbital  space  occurs 
in  the  hyperfine  interaction  H5,  given  by 

Hs  =  (2BW^  l  V1/ri  '  (36) 


where  B  is  the  Bohr  magneton,  is  the  nuclear  magneton,  is  the 
nuclear  moment,  and  I  is  the  nuclear  spin.  Now 

-  t±  -  sA  +  3r. fr. *s. 1/r?  (37) 


or 


Nq(i)  «  Aq( ± )  -  /10  J  <1(v)2(q-v)|1(q)>sv(i)C2;q_v(i)  {38) 

(we  shall  show  later  how  eq  (38)  is  obtained  from  eq  (37)). 

The  part  of  No(i)  containing  ia(i)  can  be  written  in  terms  of 
U1q,  as  in  the  second  part  of  equation  (35): 


l  V1,CJ„-v(i) 


=  <s«siisx£»c#,iu>v1  2 

2  V,q-V 


■>  r» 

A  component  of  N  =  i  can  be  written 


(39) 


Nq  =  /lU+TT  V1q  -  /To  /s  ( s+ 1 )  <£(  0)  2(  0  )  | £( 0) > 

x  l  <1(v)2(q-v)|1(q)>V1  2  . 

,,  v,q-v 


(40) 


Thus,  equation  (36)  can  be  written 


“5  -  («WI)<,/r3>  l  V? 

q 


with  Nq  given  by  equation  ( 40 ) . 
2.4  Racah  Coefficients 


(41) 


The  Racah  coefficients  arise  in  the  coupling  of  three  angular 
momenta  (Rose,  1957:  107)  to  form  a  final  resultant.  In  the  coupling  of 
the  angular  momenta,  we  consider  two  coupling  schemes: 

scheme  A:  ^  ^  -  ?12>  J12  +  l3  *  J  ,  (42) 

scheme  05  +  “  ?13'  ?13  +  ?2  “  ^  *  (43) 

Coupling  scheme  A  is  represented  by  the  wave  function 

|A>  =  l  <j1(m1)j2(m2)!  j12(m1+m2)xj12(m1+m2)j3(m3)|  j(m)> 

m1m2m3  (44) 

X  lj1m1j2m2j3m3>  7 

scheme  B  is  represented  by  the  wave  function 
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TO 


lB>  ■  I  <j1("'1)j3(m3)|  j13(m1+m3)xj13(m1+m3)j2(m2}|  j(m)> 


mim2m3  (45) 

x  I  j1m1j2m2j3m3>  . 

The  coupling  schemes  A  and  B  are  connected  by  a  unitary  transformation 

| B>  =  l  <A|B>|A>  >  (46) 

A 

the  coefficients  of  the  unitary  transformation  are  determined  by  taking 
the  inner  product  of  equation  (44)  with  equation  (45). 

We  define  the  Racah  coefficients  as  follows: 


"(hJuhih'hl)  "  [(2i12+i)(2313+i))>^  <i|B>  ‘ 


Thus, 


[(2jl2+,*)(2ji3+1)]1//2  w(323l2^13^3f^l^) 

=  I  <j1(m1)j2(m2)| j12(m1+m2)><j12(m1+m2)j3(m"m1-m2)| j(m)>  (40) 

m1m2 

x  < j1(m1 ) j3(m-m1-m2) | j13(m-m2)>< j13(m-m2J j2(m2 ) | j(m)>  . 

The  following  equation  can  be  obtained  from  equation  (48); 

<  j2(m2)Mml)*  ^12(m1+m2)><^12(m1+m2)^3(m"mrm2)l  ^(m)> 

=  J  [(2j12+l)(2j13+l)]l/2  W(j2^  1^3^12^13)  (49) 

J13 

x  <j1(m1)  j3(m-m1-m2)|  j13(m-m2)xj2(m2)  j13(m-m2)|  j(m)>  , 


which  is  a  relationship  used  often  in  our  analysis. 

The  Racah  coefficient  is  related  to  the  symmetrized  "6j"  symbol 
by  the  following  equation: 


,  .a+b+c+d  ra  b  ei 
W(abcd;ef)  =*  (-)  |d  r  f)  . 
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Certain  symmetry  relations  exist  for  the  "6j"  symbols: 


^1  j2  j3 
*“1  *2  *3 

i 


j2  j1  ^3 

z2  A1  £3 


^1  ^3  -*2 
*1  *3  *2 


j1  *2  M 

*1  j2  ^3  | 


(51) 


and  all  combinations  of  the  relations  in  equation  (47).  The  four  triads 

(j)  $2  ^3 ) '  (^1  ^2  *3^' 
a  triangle.  That  Is, 

I j1-j2l  <  j3  <  j,+j2»  (52> 

with  similar  relations  for  the  other  triads. 


( j2  £3),  and  (£^  j3)  must  be  aole  to  form 


An  example  of  the  use  of  Racah  coefficients  is  in  the  calcula¬ 
tion  of  single-electron  matrix  elements  of  the  operator 


l  <k(q-A)1(X)!k'(q)>Ck>q_A£x 


(53) 


which  arises  in  numerous  applications.  We  consider  the  matrix  element 

<rm’|E  |  £m>  =  <£(m)k’ (q)  I  £’ (m* )  X£*  IE.  ,  l£>  (54) 

K  ^  JC 

by  application  of  the  Wigner-Echart  theorem,  equation  (30).  Also,  by 
taking  the  same  matrix  element  of  equation  (53)  we  have 

<  £ '  m '  I  E  ,  |£m>  =  T  <k(q-A)1(A)|k,(q)X£'m,|C.  ,£, | £m>  .  (55) 

K  (J  K  $  CJ“  A  A 

Now  we  further  consider  the  matrix  element  in  equation  (55)  to  obtain 

<£'m'|c,  ,£.|£m>  =  J  <£'m'  |C  .  U"m,,X£"m"  |£.  I  £m>  ,  (56) 

k,q"A  A  r'm"  k,q_x  A 

where  we  have  used  matrix  algebra  on  the  product  of  two  operators.  If 
we  apply  the  Wigner-Echart  theorem  to  the  last  matrix  element  in  equa¬ 
tion  (56),  we  obtain 

<£"m"  |  £^|  £m>  *  <£(m)  1  (  A)  |  £"  (m" )  >  6U„<£  I  £1  £>  ;  (57) 

also,  m"  =  m  +  A  as  required  by  the  C-G  coefficient.  It  can  be  shown 
that 


<£«£«£>  =  /£(£+!) 


(58) 
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SE'SBiEByE'SRS  JS'WPSfflKSS?? 


ss&rnsrawssc** 


Then, 


<*"m"UAlfcn>  -  <£(m) 1 ( A)  |  A"  (m+X)  A+1 )  •  (59) 

Using  these  results  in  the  remaining  C-G  in  equation  (56),  wo  have 

crm'la  Ji(m+X)>  -  <£(m+A)k(q-A)  U’ (m*  )><&' DC  U>  .  (60) 

k,q-A  K 

Substituting  the  result  of  equations  (60)  and  (59)  into  equation  (56), 
we  have 

<rm'|C,  =  ✓*(  4+1X4'  1C.  lUXJt(m)KA)  U(m+X)> 

k,q~X  A  k 

(61) 

x  <&(m+\)k(q-\) | (m' ) >  , 

giving  the  matrix  element  in  equation  (55),  If  we  substitute  the  result 
of  equation  (61)  into  equation  (55),  then  we  have 


< Jt ' m '  |  E  ,  |*m>  =  Hi  Jt+1)  <«•'  1C.  lt>S  , 

k'  a  k 


where 


S  =  l  <k(q-A)1( A)|k‘  (q)XJt(m)  1(A)  U(m+X)X£(m+X)k(q-X)  U’  (m*  )>  .  (63) 

X 

The  last  two  C-G  coefficients  in  equation  (63)  can  be  recoupled  by  using 
equation  ( 49 )  or 

<£.(m)  1  (  A)  I  Jt(m+A)  ><fc(m+A)k(q-A)  I  (m')> 

-  I  /( 2f+1 )  ( 24+1 )  wdU'kflf  )<£(A)k(q-A)|f(q)  (g4) 


x  <£(m)f (q) | (m' )> 


The  C-G  coefficients  in  equation  (64)  can  be  rearranged  by  using  the 
symmetry  rules  to  g  i.ve 

<1( X)k(q-X) | f(q)>  =  ( - 1 ) 1 +k_f <k (q-X) 1 ( X) I f (q ) >  .  (65) 

This  C-G  coefficient  and  the  first  C-G  coefficient  in  equation  (63)  are 
the  only  two  C-G  coefficients  containing  X,  so  that 

l  <k(q-A)1(A)  |k*  (q)Xk(q-A)1(A)  |f(q)>  *  «fk,  (66) 


✓ 


because  of  the  orthogonality,  as  shown  in  equation  (19),  of  the  C-G 
coefficients.  Thus,  we  get 


-  - - 

S  *  (-1)  *  *  /(2k'+1)(2£+1)  W(£U,kiik,)<t(m)k'(q)|t,(m,)>  *  (67) 

which  when  substituted  into  equation  (62)  gives 

1 '  - 

<  fc'm'  |  E  ,  |  £m>  *  (-1)  ✓i(U1)(2)l+1)(2k'  +  1)  W(  £1£,k;£k'  )<£'  DC  ll£> 

JC  (J  K. 

X  <  £  ( m )  k '  (q)  I  £'  (in'  )>  .  (68) 


Upon  comparing  the  result  given  in  equation  (54)  with  equation  (68),  we 
have 


<£*flEklH£>  =  (-1)1+k-k' [£(£+1){2£+1)(2k,+1)  j  1/2  <£’!)Ckll£> 
x  W(  £1  £'k>  £k* )  , 


(69) 


which  is  a  useful  relation  if  we  wished  to  express  the  tensor  Ek<g  in 
terms  of  unit  spherical  tensors;  in  that  case  we  would  specialize  equa¬ 
tion  (69)  to  £'  =  £  and  simply  replace  in  equation  (34)  by  Ekg  with 
the  reduced  matrix  element  given  by  equation  (69).  We  shall  have  fre¬ 
quent  occasion  to  express  our  results  in  terms  of  Racah  coefficients  by 
using  equation  (49). 


2 . 5  Racah  Algebra 


It  is  convenient  in  many  vector  problems  to  express  the  vectors 
in  terms  of  spherical  bases  given  by 

e±1  “  *(ex  *  iey)2V2  ' 

(70) 

Then  e*  =  (-1 )ye_y  , 


e  x  e  =  -i/2<1(v)1(y)|1(u+v)>  ,  (71) 

=  6yv  * 


The  vector  A  can  be  written 


A 


A 


I  e*A 

M  U  U 


and 


=  y  e  A* 

y  U  U 


=  y  a*b 

t  u  v 
=  y  a  b* 

u  u 

y 

»  I  (-Dma  b 

44  —  U  u 

M 


Thus,  t«s  in  equation  (5)  can  be  written 


(72) 


(73) 


=  I  (74) 

so  that  the  spin-orbit  interaction  given  in  equation  (5)  is  immediately 
in  spherical  tensors,  since  and  are  spherical  tensors. 

As  an  application  of  Racah  algebra  and  some  of  the  other  mater" 
ial  discussed  above,  we  shall  derive  the  gradient  formula  (Rose,  1957: 
120).  A  convenient  form  for  the  gradient  operator  is 


and  we  would  like 


V 


s 


grad  <f>(r)Ckq(r) 


[V,*<r>ckq3 


(75) 


(76) 


First  we  observe  that 


l  <-1,X«-lCUCk , 

l  (-1)Xe  <1(0)k(0)|k'(0)xl(X)k(q)|k'(q+X)>C 

X,k*  "A 


k',q+X 


(77) 
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where  we  have  used  the  coupling  rule  for  spherical  harmonics,  equation 
(15)  (Rose,  1957:  61).  Now  we  write 

;  *  t ,  j  hi*,»  ;  .  i  .  <™> 

u6 

we  use  equation  (71)  to  eliminate  the  cross  product  to  produce 

r  x  t  m  -i/2  l  (-1)X<1(X-a)1(a)|1(X)>e,C  i  .  ,  (79) 

%  a  i "Oi  a-* A 

A/  a 

where  we  have  replaced  the  sum  on  6  by  letting  6  =  X  -  a.  Now  in  calcu¬ 
lating  the  commutation  we  need  only  consider  the  operators  in  equation 
(79);  thus,  we  need 


(80) 


[C,  i  . , C  ]  . 

L  1-a  a- A  kqJ 

Since  (J>(r)  commutes  with  £.  ,  we  need  not  consider  it  at  present. 

First  we  expand  the  commutator  to  obtain 

[ci  *  .»C.l  *  *  .C,  -  C.  C,  i  ,  ;  (81) 

1  1-a  a- A  kqJ  1-a  a-X  kq  kq  1-a  a-X 

we  then  use 

l  ”  [A  ,,C.  ]  +  C.  I  ,  (82) 

a-X  kq  1  a-X  kqJ  kq  a-X 

in  equation  (81)  to  obtain 

[VA-^Xq]  "  °1-.(Vx-Ck,l  *  Ct-aCkq*a-X  *  WA-X  '  1831 


The  last  two  terms  cancel  since  C«  and  Cv„  commute.  Thus,  we  obtain 

1-a  kq 

[C1-a*a-X'CkqJ  =  C1-at *a-X'Ckq ]  (84) 

=  C1.a/k(k+1)  <k(q)1(a-X) |k(q+a-X)>Ck  q+a_x  ,  (85) 

where  we  have  used  equation  (27)  with  La_x  =  £q_x  (which  are  identical 
in  the  commutation  brackets).  The  result  in  equation  (85)  is  not  quite 
in  the  form  we  want,  but  by  using  the  coupling  rule  for  spherical  har¬ 
monics  given  in  equation  (15),  we  finally  obtain 


[c1_aAa_x,Ckq]  =  ✓k(ktl)  <k(q)1(a-X)|k(q+a-X)>  (86) 

x  l  <k(0) 1(0) |k"(0)><k(q-a-X) 1(-a) |k"(q-X)Ck„  .  . 

k"  'q 
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In  equations  (75),  (76),  and  (79),  we  need 


(87) 


we  can  see  from  equations  (86)  and  (79)  that,  when  this  is  formed,  the 
terms  dependent  on  a  are 


S  =  l  <1  (  a-X)  1  ( -a) |l(X)><k(q)1(a-X) |k(q+a-X) > 
a 


(88) 


x  <k(q+X-a)  1  ( — ot)  |k"(q-X)>  ; 


that  is , 


[rxt,Cv  1  -  i/2  I  (-1)Ae,<k( 0)1(0) | k" ( 0 ) >/k (k+1 ) 
kq  X, a,k"  A 


(89) 


*  SCk",q-X 


with  S  given  by  equation  (88).  The  sum,  S,  given  by  equation  (88)  can 
be  reduced.  First  we  write 


<k(q) 1 ( a-X) |k(q+a-X) ><k(q+a-X) 1 (-a) |k"(q-X) > 

(90) 

«  l  /(Tf+1 ) ( 2k+iT  W(k1k" 1 ;kf )<1 ( a-X) 1 (-a) | f ( -X) ><k(q)f ( -X) | k" (q-X) >  , 

f 


where  we  have  used  equation  (49).  Thus,  the  sum  over  a  contains  the 
terms 

l  <1(a-X)1(-a) | 1(-X)><1(a-X)1(-a)|f ( -X) >  -  5f1  (91) 

a 

by  the  orthogonality  of  the  C-G  coefficients.  We  can  use  equations  (91) 
and  (90)  in  equation  (88)  to  obtain 

S  »  /3T2k+TT  W(k1k"1jk1 )<k(q) 1(-X) |k"(q-X)>  .  (92) 

Using  the  results  of  equation  (92)  in  equation  (89)  gives 

[rx$,ckq]  =  i/6k(k+1)(2k+1)  l  (-1)Ae_x  l  <k( 0 ) 1 ( 0) |k( 0) >W<k1kM 1 ;k1> 

X  k" 

x  <k(q) 1 (X) Ik"(q+X)>Ck„^q+^  ,  (93) 
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where  we  have  changed  the  sign  of  A  in  the  sum.  Multiplying  the  results 
given  in  equation  (93)  by  -i$(r)/r  and  combining  them  with  equation 
(33),  we  have  (changing  k"  to  k' ) 

,  A  f3d>  <t> 

[?.*(r)Ck  ]  =  I  (-1)*e_x  l  —  +  — /6k ( k+ I ) ( 2k-t- 1 )  w(klk’ljkl) 

A  k'L3r  r  J 

(94) 

x  <k(0)1<0)|k'(0)><k(q)1(X)!k’(q+A)>Ckl^q+x  . 


The  Racah  coefficients  in  equation  (94)  are  of  simple  form  and  are  given 
by  Rose  (1957:  227).  These  are 


W(k1k' 1 jkl ) 


_  f _ k 

|_6  ( k+ 1 ) 


(2k+1 ) 


_ K±J  .  1 l//2 

6k( 2k+1 ) 


k*  =  k  +  1 


k'  =  k  -  1 


which  are  the  only  values  of  k'  allowed.  These  results  used  in  equation 
(94)  can  be  written  as 

[v,<fr(r)Ckq]  =  l  (-1 ) Ae_x  |t<k(q)1(X)|k*(q+X)>Ck,fq+xDk’(Mr)  ,  (96) 

where 


k'  ^  /  k+1  /3_  k\ 

V  2k+ 1  \3r  ~  r/  ' 

D*'  .  .  +  *±±\ 

D  \2k+1  \3r  r  /  ' 


k’  =  k  +  1  , 


k1  =  k  -  1  , 


and  we  have  used  the  result 


<k(0)1(0)|k+1(0)>  =  and 


<k(0) 1(0) | k— 1 ( 0 ) >  =  - 


(96a) 


(96b) 


(from  Rose,  1957:  225).  The  two  most  common  forms  of  d>(r)  that  we  will 

If  ]/X  1  If 

encounter  are  r  and  1/r  .  For  <f>(r)  =  r  ,  we  obtain 


[Vy,rkCkq]  =  -  /k(2k+1)  <k(q)1(u)|k-1(q+U)>rk'1Ck_^q+ij  ,  (97) 
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and  for  $  =*  1/r 


k+1 

/ 


K 


,  /  ^+1 
1/r  C 


-  /(k+1 ) ( 2k+ 1 )  <k(q)1(u) |k+1{q+M)>1/rk+2C 


kq  ■ 


k+1 ,q+u 


(98) 


The  results  given  in  equation  (98)  are  easily  checked  for  k  =  0,  since 
for  k  *  0  we  have 


[V  /  1/r  ]  =  -<0(q)1(M)H(q+y)>1/r2Ci  ,  (98a) 

W  i  #q+M 


and  from  the  properties  of  the  C-G  coefficients,  we  know  that  q  =*  o 
and  <0(p)1(p)|1('rj)>  =  1.  Then 

[VVr'C1q]  =  *  ("b) 

Also,  we  know  from  vector  analysis  that 

grad  1/r  =  -r/r3  (98c) 


and 


#  -  r  l  (-D^C^  .  (98d) 

Then  we  substitute  equation  (98d)  in  (98c)  to  obtain 

(grad  1/r)^  =  -C^/r2  ,  (98e) 

which  is  identical  with  the  result  of  equation  (98b).  We  shall  use  the 
result  given  in  equation  (94)  frequently  later  on,  particularly  in  the 
form  given  in  equations  (97)  and  (98). 


3.  FREE- ION  HAMILTONIAN 

3.1  Background  for  Free  Rare-Earth  Ions 


The  approximations  made  in  the  analysis  of  rare-earth  ions  are 
not  new.  In  fact,  they  go  back  to  the  old  Bohr  orbit  theory.  Since 
many  of  you  may  not  be  familiar  with  these  assumptions  and  may  not 
remember  many  of  the  concepts  and  much  of  the  technical  jargon  used  in 
the  field  of  atomic  spectra,  I  will  review  some  of  these  briefly. 
Although  generalities  may  exist,  I  will  stick  strictly  to  those  concepts 
which  apply  to  rare-earth  ions  or,  more  strictly,  triply  ionized  rare- 
earth  ions.  The  triply  ionized  rare-earth  ions  are  characterized  by  the 
electronic  structure  shown  in  table  1 . 
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TABLE  1.  ELECTRONIC  STRUCTURE  OF  TRIPLY  IONIZED 
RARE-EARTH  IONS 


Number  Element- 

Symbol 

Outermost  electron  shell 

57 

Lanthanum 

La 

4d104f °5s25p6 

58 

Cerium 

Ce 

4d104f l5s25p6 

59 

Praseodymium 

Pr 

4dl04f25s25p6 

60 

Neodymium 

Nd 

4d104f 3Ss2$p6 

61 

Promethium 

Pm 

4d104f 45s25p6 

62 

Samarium 

Sm 

4d104f55s25p6 

63 

Europium 

Eu 

4d104f65s25p6 

64 

Gadolinium 

Gd 

4d104f75s25p6 

65 

Terbium 

Th 

4d104f85s25p6 

66 

Dysprosium 

Dy 

4d104f95s2Sp6 

67 

Holmium 

He 

4d104f 105s25p6 

68 

Erbium 

Er 

4d' °4f 1 15s25pb 

69 

Thulium 

Tm 

4d104f 125s25p6 

70 

Ytterbium 

Yb 

4d104f 135g25p6 

71 

Lutetium 

Lu 

4d104f 145s25p6 

In  the  rare-earth  series ,  it  is  assumed  that  the  atomic  inter¬ 
actions  are  very  strong;  thus,  when  an  ion  is  placed  in  a  crystal,  the 
crystalline  electric  field  acts  as  a  perturbation  on  the  ion.  This 
assumption  allows  the  notation  developed  for  the  free  ion  to  be  used 
with  the  reservation  that  many  of  the  "good"  quantum  numbers  of  the  free 
ion  are  not  quite  good  when  the  ion  is  present  in  the  crystal.  It  is 
assumed  that  the  free  rare-earth  ions  have  the  zeroth-order  Hamiltonian 

N 


2m 


(99) 


where  p^  is  the  momentum  of  the  ith  electron  and  Ufr^)  is  an  appropriate 
spherical  average  potential  of  the  remaining  (other  than  the  N,  4f) 
electrons  in  the  ion.  The  single-electron  solutions  to  equation  (99) 
are  taken  in  the  form  (Schiff,  1968) 

*  ■  R„jlr”Wr'>  • 

A 

where  the  Y^m(r)  are  the  spherical  harmonics  with  £  -  3  for  f  electrons. 
(Remember  that  £  «  0  for  s,  £  -  1  for  p,  and  £  =  2  for  d  electrons.) 
The  radial  functions  in  equation  (100)  are  taken  to  be  the  same  for  all 
the  f  electrons  in  the  ion,  while  the  angular  functions,  along  with  the 
spin  of  each  electron,  must  form  a  determinantal  function  so  as  to  obey 
the  exclusion  principle.  Depending  on  whatever  determinantal  function 
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is  chosen,  the  radial  functions  can  be  found  by  some  self-consistent 
method.  These  radial  functions  (Freeman  and  Watson,  1962 j  Fraga  et  al, 
1976?  Cowan  and  Griffin,  1976)  have  been  found  for  the  Hund  ground 
states  of  all  the  rare  earths  from  cerium  (Ce3+)  through  ytterbium 
(Yb3+). 


The  Hund  ground  state  for  the  rare-earth  ions  with  N  <  7  is 
given  by  assuming  that  all  N  spins  are  parallel  and  that  each  angular 
momentum  projection  is  the  maximum  allowed  by  the  exclusion  principle. 
(In  eq  (100),  l  is  the  angular  momentum  and  m  is  its  projection.)  Thus, 
the  Hund  ground  state  for  two  electrons  is  the  determinantal  function 

«(1)*£,|(1)«(2)Yti£.1(2)  -  a(1)Y^£_1(1)a(2)Y£^(2)  , 


where  a  is  the  spin  "up"  wave  function  (B  =  spin  "down").  A  convenient 
notation  for  such  a  determinant  is 

{£,£-1}  ,  (101) 

where  the  upper  sign  is  the  spin  projection  (+  =  up,  -  =  down)  and  l  and 
£-1  are  the  z  projection  of  the  angular  momentum  (m  in  eq  (100)).  Thus, 
the  Hund  ground  state  equation  (101)  for  4f^  praseodymium  (Pr3+)  has 
total  spin,  S,  and  total  angular  momentum,  L,  given  by 

S  =  1/2  +  1/2  «  1, 

L=£+£-1=2£-1=5  (for  f,  i  -  3)  . 

Hence,  the  ground  state  is  L  =  5,  with  multiplicity  of  2S  +  1  =  3.  In 
the  so-called  Russell-Saunders  notation,  this  state  is  referred  to  as 
3H,  as  given  by  the  table: 

Total  angular  momentum,  L,  of  ion  01234567.... 

Russell-Saunders  label  for  the  state  SPDFGHIK.  .  .  . 

(continues  alphabetically). 

In  this  notation,  the  technical  reference  to  such  a  state  is  a  term; 
other  terms  for  Pr3  +  are  ^I,  3P,  and  (Condon  and  Shortly,  1959).  For 
the  ion  Ce3+,  which  has  one  f  electron,  the  atom  notation  becomes  iden¬ 
tical  to  that  of  the  ion,  that  is,  i  =  L  =  3  and  s  ■  S  =  1/2  with  the 
2 

single  term  F.  Those  ions  in  the  series  for  N  >  7  have  the  same  terms 
as  for  N  <  7,  and  their  Hund  state  can  be  constructed  simply  as 

{£,$-1,1-2,  .  .  . ,£-(2*+1);£,Z-1,  .  .  .  ,Z-p+1}  ,  (102) 

N 

where  the  number  of  f  electrons  is  N  =  7  +  p.  The  i  shell  is  com¬ 
pletely  filled  when  N  =  2(2£+1),  which  for  f  electrons  is  triply  ionized 
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lutetium  (Lu^+).  As  an  example,  consider  triply  ionized  terbium  (Tb^+), 

Q 

which  has  the  4f  configuration.  The  determinantal  wave  function  is 

t  f  +•  +  +  +  +  - 

{  3  2  1  0  -1  -2  -3  3}  ,  (  103) 

where  total  angular  momentum  L  =  3  and  total  spin  angular  momentum 
S  »  6/2  =3.  Thus,  the  Hund  ground  state  is  'F.  In  all  cases,  the  Hund 
term  has  been  found  to  have  the  lowest  energy  in  atomic  systems.  In 
general,  the  wave  functions  for  the  higher  terms  are  very  difficult  to 
construct,  but  sophisticated  techniques  have  been  devised  for  the 
orderly  development  of  a  sat  of  wave  functions  for  each  ion  in  the 
entire  rare-earth  series  (Nielson  and  Koster,  1964).  The  Hamiltonian 
given  in  equation  (99)  has  the  same  value  for  all  terms  of  the  configu¬ 
ration  4fN;  consequently,  we  ignore  HQ  in  the  future  discussion. 

3 . 2  Significant  Free  Rare-Earth  Ion  Interactions 
3.2.1  Coulomb  Interaction 


The  largest  contribution  to  the  Hamiltonian  for  a  rare-earth 
ion  is  the  electrostatic  interaction  of  the  4f  electrons,  which  may  be 
written 


H 


1 


(  104) 


where 


The  matrix  elements 
term)  of  Pr3+  are 


*ij  =  ?i"*j  ' 

of  this  interaction  for  the  state  {5  h  (the 


<3H|H.|3H>  Fn  -  25F  -  5  IF.  -  1 3F 

1  0  2  4  6 


(  105) 


where  the  F^  are  frequently  referred  to  as  the  Slater  parameters.  These 
F^  are  related  to  the  parameters  F^k-  (also  frequently  referred  to  as 
"Slater  integrals")  by  the  following; 

F  -  f(0),  =  F<2)/225,  Fa  »  Ft4)/1089,  and  Fc  =  25F(6)/184,041  . 

0  2  4  b 


(This  relation  is  strictly  for  f  electrons;  a  different  relation  exists 
for  d  electrons.)  The  F'*'  are  radial  expectation  values  given  by 


,(k) 


eZ  /o  Q  “T+T  fRn?(rl)RnJr2^2  drl  dr2 

r> 


(  106) 


/ 
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where 


Q  R^(r)  ir  -  ,  , 


r  .  and 
3 


The  confusion  created  by  the  two  sets  of  F's  does  not  end  at  this  point; 
there  is  still  one  further  set  of  parameters,  Ek,  called  the  Racah 
parameters,  which  are  used  to  express  the  matrix  elements  of  equation 
(104).  The  relation  of  Ek  to  F^k'  and  F^  is  given  in  a  number  of  places 
(for  example,  Judd,  1963:  206). 

3.2.2  Spin-Orbit  Interaction 


The  second  interaction  of  reasonable  magnitude  in  the  free 
ion  is  the  spin-orbit  coupling,  which  is 


r< 

r> 


—  if  r.  < 

r  j  i 


-t  if  r.  > 
ri  1 


N 

H-  =  l  5(r  )I .«s  ,  (107) 

i=1  1  1  1 

where 

«(*i) 

This  interaction  was  derived  from  relativity  theory  in  the  Bohr  orbit 
quantum  mechanics,  but  is  a  natural  consequence  of^  a  nonrelativistic 
approximation  to  the  Dirac  equation.  Values  of  F  '  and  x,  (where  x,  = 
<4f | £(r) | 4f > )  from  Hartree-Fock  wave  functions  are  given  in  tables  2  and 
3.  In  the  rare-earth  series,  the  interaction,  H 2,  is  quite  strong  and 
is  in  general  much  larger  than  the  interaction  of  the  rare-earth 
electrons  with  the  crystal  fields. 

Consequently,  it  is  convenient  to  perform  all  the  calcu¬ 
lations  in  a  set  of  basis  functions  in  which  H1  and  H2  are  diagonal. 
The  set  of  functions  that  achieves  this  is  the  total  angular  momentum 
function  |JMj>,  where  the  total  angular  momentum  operator  J  55  L  + 


du(r.) 


„  2  2  r, 
2m  c  1 


dr , 
x 
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TABLE  2.  RELATIVISTIC  HARTREE-POCK  INTEGRALS  FDR  TRIPLY 
IONIZED  RARE-EARTH  IONSa  (all  in  cm"1) 


N 

r3+ 

F2 

F4 

F6 

C 

M° 

M2 

M4 

1 

Ce 

_ 

696.41 

_ 

2 

Pr 

98723 

61937 

44564 

820.22 

1.991 

1.110 

0.752 

3 

Nd 

102720 

64462 

46386 

950.51 

2.237 

1.248 

0.846 

4 

Pm 

106520 

66856 

48111 

1091.46 

2.492 

1.391 

0.943 

5 

Sm 

110157 

69143 

49758 

1243.60 

2.756 

1.540 

1.044 

6 

Eu 

113663 

71373 

51342 

1407.71 

3.G31 

1.694 

1.149 

7 

Gd 

117058 

73470 

52873 

1584.45 

3.318 

1.855 

1.258 

8 

Tb 

120366 

75541 

54361 

1774.46 

3.615 

2.022 

1.372 

9 

Dy 

123592 

77558 

55810 

1998.44 

3.924 

2.195 

1.490 

10 

Ho 

126751 

79530 

57227 

2197.06 

4.246 

2.376 

1.612 

11 

Er 

129850 

81462 

58615 

2431.00 

4.580 

2.563 

1.739 

12 

Tm 

132897 

83361 

59978 

2680.97 

4.928 

2.758 

1.872 

13 

Yb 

“ 

2947.69 

- 

aR. 

D,  Cowan  and  D.  C, 

Griffin, 

1976,  J. 

Opt.  Soc. 

Am.  66/  10 10 • 

TABLE  3.  NONRELATIVISTIC  HARTREE-FOCK  INTEGRALS  FOR  TRIPLY  IONIZED 
RARE-EARTH  IONS3  (all  in  Cm-1 ) 


N 

r3+ 

F2 

F4 

F6 

C 

M° 

M2 

M4 

1 

Ce 

. 

_ 

778.14 

2 

Pr 

1051 20 

66213 

47718 

919. 16 

2.2610 

1 . 2669 

0.8601 

3 

Nd 

109731 

69165 

49860 

1069.87 

2.5478 

1.4294 

0.9709 

4 

Pm 

1 1 3640 

71641 

51647 

1228.24 

2.8261 

1.5865 

1.0779 

5 

Sm 

117222 

73893 

53269 

1397.79 

3.1079 

1.7453 

1.1859 

6 

Eu 

120885 

76204 

54937 

1583.54 

3.4093 

1  .9153 

1 .3016 

7 

Gd 

124644 

78585 

56655 

1786.68 

3.7320 

2.0976 

1.4258 

8 

Tb 

127137 

80091 

57722 

1990.51 

4.01  30 

2.2544 

1 .5320 

9 

Dy 

129960 

81829 

58962 

2214.87 

4.3231 

2.4281 

1 .6500 

10 

Ho 

132929 

83670 

60281 

2458.58 

4.6552 

2.6144 

1.7766 

11 

Er 

135859 

85486 

61580 

2719.76 

5.0004 

2.8081 

1.9081 

12 

Tm 

138754 

87276 

62864 

2999.22 

5.3590 

3.0092 

2.0448 

13 

Yb 

- 

- 

3299.82 

aS.  Fraga,  K.  M .  S.  Saxena,  and  J.  Karwowski,  1976,  Physical  Sciences 
Data  5/  Handbook  of  Atomic  Data  (Elsevier ,  New  York), 
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The  spin-orbit  interaction  Hj  given  in  equation  (107)  com¬ 
mutes  with  the  total  angular  momentum  and,  consequently,  since  H.,  also 
commutes  with  J2,  the  wave  functions  can  be  characterized  by  the  eigen¬ 
values  of  J2  and  J_.  That  is,  we  can  write  or  (JM>  for  the  wave 

z  JM 

functions  with 

J2|JM>  ■  J(J  +  1 ) | JM>  (108) 


and 


JZ|JM>  =  M|JM>  . 

For  any  term  of  given  L  and  S  (eigenvalues  of  L2  and  §2),  the  values  of 
J  are  restricted  to 


|L  -  S|  <  J  <  |L  +  S|  . 


Then  the  wave  functions  are  customarily  written 
have 


*JMLS  °r  ,JMLS>' 


and  we 


L2|JMLS>  «  L(L  +  1 ) f JMLS>  , 
S2  j JMLS>  =  S(S  +  1 ) | JMLS>  , 
J2|JMLS>  =  J(J  +  1 )  j  JMLS>  , 
JZ|JMLS>  =  M | JMLS>  , 


(109) 


and 


CJ’M'L'S’ | H ^  +  H2|JMLS>  «  •  (HO) 


As  implied  In  equation  (110),  the  energy  and  Hj  is  independent  of  M, 
or  each  J  level  of  the  free  ion  is  2J+1-fold  degenerate.  The  matrix 
elements  in  equation  (110)  do  not  vanish  generally  for  L'  =*  L  ±  1  and  S' 
■  S  ±  1j  thus,  L  and  S  are  not  strictly  good  quantum  numbers.  Neverthe¬ 
less,  the  energy  levels  are  labeled  as  though  they  were,  as  in  the 
Russell-Saunders  notation,  2S+^Lj.  An  example  of  the  energy  levels  for 
the  4f2  configuration  of  the  free  ion  (Pr3+)  is  given  in  table  4.  Also 
included  is  the  same  ion,  Pr3+,  in  the  host  materials  lanthanum  tri¬ 
chloride  (LaCl^)  and  lanthanum  trifluoride  (LaF^). 

The  results  in  table  4  are  interesting  in  that  they  show  that 
most  of  the  energy  levels  observed  in  the  free  ion  are  lowered  when  the 
ion  is  embedded  in  a  solid.  This  shift  In  the  energy  levels  is  a 
general  effect  and  is  not  restricted  to  Pr3+,  but  exists  in  all  the 
rare-earth  ions  where  a  comparison  with  energy  levels  of  the  free  ion 
can  be  made.  In  fact,  this  shift  has  been  observed  by  Low  (  1958a, b)  in 
ions  with  an  unfilled  d  shell.  The  first  explanation  of  this  shift  in 
energies  was  by  Morrison  et  al  (1967),  where  it  was  shown  that  if  the 
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ion  under  investigation  was  assumed  to  be  embedded  in  a  solid  of  homo¬ 
geneous  dielectric  constant,  e,  then  a  decrease  in  the  Slater  integrals 
is  given  by 

AF(k)  =  -e2(e-1)(k+1)(<rk>)2/{b2k+1[e+Me+1)]  }  ,  (111) 

where  b  is  the  radius  of  a  fictitious  cavity  surrounding  the  rare-earth 
ion.  The  result  given  in  equation  (111)  was  first  successfully  applied 
to  Co++  in  MgAl204.  Later,  Newman  (1973)  showed  that  the  shift  in 
given  in  equation  (111)  was  sufficiently  large  to  predict  the  shifts  in 
the  energy  levels  for  rare-earth  ions.  More  recently,  Morrison  (1980) 
derived  the  result 


AF 


(k) 


Vi6 

2k+4 


2 

9 


(112) 


k  • 

where  is  the  polarizability  of  the  Z^  ligands  at  and  <r  >  is  the 
radial  expectation  value  of  r  .  The  result  given  in  equation  (112)  is 
believed  to  be  more  fundamental  than  that  of  equation  (111)  because  the 
latter  explicitly  accounts  for  the  local  coordination  of  the  rare-earth 
ion.  Morrison  (1980)  gives  a  predicted  shift  in  the  spin-orbit  param¬ 
eter,  £,  but  because  of  the  smallness  of  the  predicted  shift  and  the 
errors  in  the  fitting  of  the  experimental  data,  no  comparisons  were 
made. 


TABLE  4.  FREE  ION  ENERGY  LEVELS  OF  TRIPLY  IONIZED  PRASEODYMIUM 
AND  CORRESPONDING  CENTROIDS  IN  TWO  CRYSTALS3  (all  in  cm-1) 


[LS]J 


3 

3 

3 

3 

3 

3 

1 

1 

3 

3 

1 

3 

1 


Free 

LaCl3 

La  Fj 

0 

0 

0 

2152 

2119 

2163 

4389 

4307 

4287 

4997 

4848 

5015 

6415 

6248 

6368 

6855 

6684 

6831 

9921 

9704 

9801 

17334 

16640 

16847 

21390 

20385 

20727 

22007 

21987 

21314 

22212 

21327 

- 

23160 

22142 

22546 

50090 

48710 

46786 

aG.  H.  Dieke, 
Ions  in  Crystals 


1968,  Spectra  and  Energy  levels  of  Rare  Earth 
( Interscience ,  New  York),  p.  200. 
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Because  of  the  lack  of  experimental  data  on  the  free-ion 
spectra  of  rare-earth  ions,  measurement  of  the  shift  in  the  Slater 
Integrals  is  possible  only  for  Pr  +  .  The  experimental  F'  for  triply 
ionized  rare  earths  in  LaF^  and  LaCl^  have  been  obtained  by  Carnall  et 
al  (  1978),  and  these  results  are  given  in  tables  5  and  6.  These  data 
can  be  used  in  conjunction  with  equation  (112)  to  obtain  results  that 
can  perhaps  be  applied  to  an  arbitrary  host  material  to  predict  a  priori 
the  energy  level  shift  of  that  host. 

TABLE  5.  FREE-ION  PARAMETERS  FOR  TRIPLY  IONIZED  RARE-EARTH  IONS 
IN  LaClj  OBTAINED  FROM  FITTING  EXPERIMENTAL  DATA3  (all  in  cm-1) 


Ionb 

F2 

F4 

F6 

a 

6 

Y 

c 

Pr 

68368 

50008 

32743 

22.9 

-674 

(  1520) 

744 

Nd 

71866 

52132 

35473 

22.1 

-650 

1586 

880 

Pm 

75808 

54348 

38824 

21.0 

-645 

1425 

1022 

Sm 

78125 

56809 

40091 

21.6 

-724 

(  1700) 

1 168 

Eu 

84399 

60343 

41600 

16.8 

(-640) 

(  1750) 

1331 

Gd 

85200 

60399 

44847 

(  19) 

(-643) 

1644 

(  1513 

Tb 

90012 

64327 

42951 

17.5 

(-630) 

(  1880) 

1707 

Dy 

92750 

65699 

45549 

17.2 

-622 

1881 

1920 

Ho 

95466 

67238 

46724 

17.2 

-621 

2092 

2137 

Er 

98203 

69647 

49087 

15.9 

-632 

(2017) 

2370 

aW.  T.  Carnall,  H.  Crosswhite,  and  H.  M.  Crosswhite ,  1973,  Argonne 


National  Laboratory,  ANL-78-XX-95 . 

Values  in  parentheses  were  not  varied  in  the  fitting. 


TABLE  6.  FREE-ION  PARAMETERS  FOR  TRIPLY  IONIZED  RARE-EARTH  IONS  IN 

r  nDTlATUPn  CDAU  PTITHHTHC  PVn.OTU.UTIT  1  1  4n  ^  \ 


LaF3 

OBTAINED  FROM 

FITTING 

EXPERIMENTAL  DATA3  (all  in  cm-1) 

Ionb 

F2 

F4 

F6 

a  8  Y 

C 

Pr 

69305 

50675 

32813 

(21) 

-842 

1625 

750.8 

Nd 

73036 

52624 

35793 

21.28 

-583 

1443 

884.9 

Pm 

(77000) 

(55000) 

(37500) 

(21.00) 

(-560) 

(  1400) 

(1022.) 

Sm 

79914 

57256 

40424 

20.07 

-563 

1436 

1177.2 

Eu 

(84000) 

(60000) 

(42500) 

(20.0) 

(-570) 

(1450) 

(1327) 

Gd 

85587 

61361 

45055 

(20.0) 

(-590) 

(1450) 

1503.5 

Tb 

91220 

65798 

43661 

19.81 

(-600) 

(1400) 

1702 

Dy 

94877 

67470 

45745 

17.64 

-608 

1498 

1912 

Ho 

97025 

68885 

47744 

18.98 

-579 

1570 

2144 

Er 

100274 

70555 

49900 

17.88 

-599 

1719 

2381 

Tm 

102459 

72424 

51380 

(17.0) 

-737 

(1700) 

2640 

aW.  T .  Carnall ,  H.  Crosswhite,  and  H.  N .  Crosswhite,  1978,  Argonne 


National  Laboratory ,  ANL-78-XX-9 5. 

b Values  in  parentheses  were  not  varied  in  the  fitting . 
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3.2.3  Interconfigurational  Interaction 

An  interaction  that  has  been  frequently  used  in  fitting  the 
"free"  ion  levels  of  rare-earth  ions  in  a  crystal  is  the  so-called 
interconfigurational  mixing  or  the  Trees  interaction.  This  interaction 
has  been  parametrized  by  Wybourne  and  Rajnak  (Wybourne,  1965)  and  is 

H1()  =  oL(L+1)  +  0G(G2)  +  YG(r?)  ,  (113) 

where  a,  0,  and  y  are  parameters  adjusted  to  fit  the  experimental  data. 
The  operator  G(G2)  is  the  Casimir  operator  for  the  group  G2,  and  G (Gy) 
is  the  similar  operator  for  R7  (note  that  L2  *  L(L+1)  is  the  Casimir* s 
operator  for  the  group  R3).  The  values  for  these  operators  for  all  the 
states  are  tabulated  by  Wybourne  (1965:  73).  The  values  for  the  state 
of  f2  are  given  in  table  7.  The  values  of  a,  8,  and  y  obtained  by 
fitting  experimental  data  for  the  rare-earth  ions  are  given  in  tables  5 
and  6.  To  my  knowledge,  no  successful  attempts  to  derive  theoretical 
values  of  o,  8,  and  y  have  been  published. 


3.2.4  Other  Interactions 

Many  other  interactions 
are  considered  in  the  free  ion,  such 
as  spin-other-orbit,  orbit-orbit, 
and  configuration  interaction.  All 
of  these  to  a  greater  or  lesser 
extent  improve  the  fit  of  theoret¬ 
ical  energy  levels  to  the  experi¬ 
mental  data.  We  will  omit  these 
interactions  from  further  discussion 
since  H1 ,  H2,  and  H1Q  give  a 

sufficient  representation  of  the 
free  ion  for  our  purposes  here. 
However,  we  shall  list  a  number  of 
interactions  including  the  above 
which  have  been  considered  by 
various  research  workers  (Wortman  et 
al,  1973a, b): 

H1  =  the  Coulomb  interaction 
H2  =  the  spin-orbit  interaction 
Hj  =  the  crystal-field  interaction 
H4  =  the  interaction  with  a  magnetic 
field  (Zeeman  effect) 

H5  *  the  hyperfine  interaction 
Hg  *  the  spin-spin  interaction 

The  notation  listed  above  is  that  of 
additions . 


TABLE  7.  EIGENVALUES  OF 
CASIMIR* S  OPERATORS  FOR 
STATES  OF  f2 


State 

a 

12b 

5y 

3P 

2 

12 

5 

3F 

12 

6 

5 

3h 

30 

12 

5 

’s 

0 

0 

0 

’d 

6 

14 

7 

’g 

20 

14 

7 

’i 

42 

14 

7 

Hy  =  the  nuclear  quadrupole 
interaction 

Hg  =  the  spin-other-orbit 
interaction 

Hg  =>  the  orbit-orbit  interaction 
H1 0  *  the  interconfigurational 
interaction 

H1 1  -  the  spin-crystal-field 
interaction 

Judd  (1963),  with  a  few  obvious 
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3.3  Summary 


We  have  considered  the  Coulomb  interaction,  H^,  and  the  spin- 
orbit  interaction,  H2,  for  the  configuration  4fN  in  the  free  ion.  The 
wave  functions  that  are  chosen  as  a  basis  for  diagonalization  of  these 
interactions  is  |JMLS>,  and  the  resulting  energy  levels  are  labeled 
according  to  the  Russel 1-Saunders  notation  as  given  in  section  3.1. 
This  same  notation  (plus  additional  quantum  numbers)  will  be  used  for 
describing  an  ion  in  a  crystal.  The  values  of  <r^>  that  are  needed  in 
equations  (111)  and  (112)  are  given  in  table  8.  The  wave  functions  used 
for  the  calculation  of  the  energy  levels  of  a  rare-earth  ion  in  a  solid 
will  be  the  combination  that  simultaneously  diagonalizes  H1  and  H2« 


TABLE  8.  NONRELATIVISTIC  HARTREE-FOCK  EXPECTATION 
VALUES  OF  <rk>,  NUCLEAR  SPIN,  AND  MAGNETIC  MOMENTS  OF  TRIPLY 
IONIZED  RARE-EARTH  IONS3  (atomic  units) 


N 

r3+ 

<r2> 

<r4> 

A 

<J\ 

V 

<r“3> 

I 

1 

Ce 

1.1722 

3.0818 

15.549 

4.88571 

7/2 

0.9 

2 

Pr 

1.0632 

2.5217 

11.492 

5.52708 

5/2 

4.3 

3 

Nd 

0,97822 

2.  1317 

8.9525 

6. 17823 

7/2 

-1.08 

4 

Pm 

0.91401 

1.8701 

7.4224 

6.82779 

7/2 

2.7 

5 

Sm 

0.86059 

1.6700 

6.3365 

7.49427 

7/2 

-0.67 

6 

Eu 

0.81064 

1.4898 

5.3886 

8.20108 

5/2 

3.464 

7 

Gd 

0.76368 

1.3267 

4.5589 

8.95054 

3/2 

-0.254 

8 

Tb 

0.73523 

1.2508 

4.2673 

9.66059 

3/2 

1.99 

9 

Dy 

0.70484 

1.1644 

3.9002 

10.42528 

5/2 

-0.46 

10 

Ho 

0.67481 

1.0788 

3.5296 

11.23504 

7/2 

4.12 

11 

Er 

0.64714 

1.0028 

3.2111 

12.07835 

7/2 

-0.564 

12 

Tm 

0.62154 

0.93514 

2.9355 

12.95607 

2 

0.047 

13 

Yb 

0.59678 

0.87030 

2.6713 

13.87743 

5/2 

-0.6776 

aS. 

Fraga,  K.  H. 

S.  Saxena, 

and  J.  Karwowski,  1 976, 

Physical 

Sciences 

Data  5, 

Handbook  of 

Atomic  Data 

( Elsevier  t 

New  York), 

4 .  CRYSTAL-FIELD  INTERACTIONS 

4. 1  Phenomenological  Theory  of  Crystal  Fields 

In  the  presence  of  a  crystal  field,  we  take  the  interaction  of 
the  rare-earth  ion  whose  configuration  is  4fN  as 

N 

J,  <WJi)  '  <"4> 
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H3  ”  l  Bkq 
kq 


V 


where  the  are  unnprmalized  spherical  harmonics  given  by 

Ckq(r)  =  Ykq(r)  * 

The  use  of  in  place  of  Y.  seems  trivial,  but  it  saves  writing  a 
tremendous  number  of  factors  of  4tt  and  2k+1  in  the  interactions.  The 
use  of  the  in  expressions  for  electronic  interactions  (along  with 
other  shorthand  notation  that  we  will  not  use)  is  practically  universal. 
The  number  of  terms  in  equation  (114)  that  needs  to  be  considered  is 
limited  by  the  symmetry  of  the  site  occupied  by  the  rare-earth  ion. 
Also,  since  we  will  be  discussing  only  the  4fN  configuration,  k  is 
limited  to  values  of  6  or  less.  This  limitation  arises  because,  inde¬ 
pendent  of  the  basis  chosen,  individual  matrix  elements  of  Ct  will  have 
to  be  considered,  and  these  are  such  that  =  0  if  k  >  2f  =  6. 

The  highest  possible  symmetry  that  an  ion  can  experience  in  a 
crystal  is  full  cubic  (for  example,  Cs2NaLnClg,  where  Ln  =  any  rare- 
earth  ion);  in  that  case,  two  are  necessary,  B4Q  and  B^q,  since  B44 
=  (5//7o)B4p  and  B&4  =  -  /( 7/2 )  and  =  0  for  k  odd.  Because  of 
the  simplicity  of  the  calculations,  the  ion  is  often  assumed  to  occupy  a 
point  of  cubic  symmetry  with  perturbations.  In  general,  this  assumption 
is  not  very  productive.  On  the  other  hand,  in  many  solids  the  rare- 
earth  ions  occupy  a  site  with  no  symmetry  (such  as  LnP5C>14,  where  Ln  =  a 
rare-earth  ion;  these  are  the  pentaphosphates ,  one  of  the  best  known 
laser  host  materials  for  Nd).  In  this  case,  all  the  B,  are  allowed, 
which  for  the  4fN  configuration  are  B2Q,  B21,  B22,  B4Q,  b41,  .  .  . 
b66*  Since  each  ,  q  *  0,  is  complex,  there  are  27  parameters. 
Unfortunately,  the  technical  importance  of  the  rare  earths  in  tnese  low- 
symmetry  solids  complicates  the  problem  of  calculating  energy  levels  for 
practical  application.  At  present,  we  assume  that  the  ion  is  in  a 
slightly  more  symmetric  position  in  the  pentaphosphates  for  simplicity 
in  the  calculation  (Morrison,  Wortman,  Karayianis,  1977)  but  have  no 
idea  how  good  these  approximations  are . 

We  have  been  assuming  in  equation  (114)  that  the  B^  are  the 
same  for  each  electron,  which  is  a  universal  assumption  made  when  deal¬ 
ing  with  rare-earth  ions.  We  will  continue  with  this  assumption  here. 
The  are  referred  to  as  the  crystal-field  parameters  and,  as  the  name 
implies,  are  used  as  parameters  when  fitting  the  experimental  data. 

In  the  early  days  of  the  experimental  investigations  of  rare- 
earth  ions  in  solids,  relatively  few  materials  could  be  grown  as  single 
crystals  with  a  small  amount  of  rare  earth  contained  substitutionally . 
An  exception  was  lanthanum  trichloride  (LaCl^),  which  was  the  crystal 
used  by  the  Johns  Hopkins  Laboratory  under  the  direction  of  Gerhardt 
Diexe,  who  in  the  span  of  10  to  15  years,  along  with  his  graduate  stu¬ 
dents,  had  reported  the  absorption  and  emission  spectra  of  nearly  every 
rare-earth  ion.  In  one  of  Dieke's  publications  (Dieke  and  Crosswhite, 
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1963),  he  gave  an  energy  level  diagram  of  the  results  of  this  research. 
This  diagram  became  so  popular  that  it  is  generally  referred  to  as  a 
"Dieke  chart"  and  is  found  in  practically  every  laboratory  where  the 
spectra  of  rare-earth  ions  are  being  investigated  or  rare-earth  ion 
lasers  are  being  designed.  These  Dieke  charts  are  invaluable  for  the 

investigation  of  a  new  crystalline  material  (commonly  referred  to  as  a 

"host")  that  has  been  grown  to  contain  a  low  atomic  percentage  of  a 
particular  rare-earth  ion.  A  common  method  for  referring  to  a  substi¬ 
tutional  doping  is  *3  ( i_x )Nd3xA^5012 *  for  x  x  100  percent  neodymium- 

doped  yttrium  aluminum  garnet.  When  a  solid  contains  a  rare-earth  ion 
as  one  of  its  constituents,  such  as  the  excellent  laser  material 
NdPgO^,  it  is  referred  to  as  a  stoichiometric  laser  material.  A  pri¬ 
mary  requirement  of  a  host  material,  besides  being  transparent  in  the 
region  of  interest,  is  that  cne  of  the  atomic  constituents  be  trivalent 
and  have  nearly  the  same  ionic  radius  as  the  ion  in  the  rare-earth 

series.  Also,  the  solid  should  have  no  absorption  bands  in  the  region 
of  interest.  The  most  common  ion  meeting  these  requirements  is  yttrium, 
and  most  of  the  rare-earth  doped  laser  crystals  have  this  element  as  a 

constitutent,  although  several  of  the  rare-earth  ions  themselves  meet 

3+  3+ 

these  requirements,  such  as  La  ,  Lu  ,  and  for  particular  purposes 
gadolinium  Gd3+. 

4.1.1  Matrix  Elements  of  H3  in  J  States 


In  order  to  make  full  use  of  tabulated  data  in  our  calcula¬ 
tions,  it  is  necessary  to  make  some  modifications  in  equation  (114). 
Neilson  and  Koster  (1964)  have  calculated  the  reduced  matrix  elements  of 
the  unit  spherical  tensors  introduced  by  Racah.  The  Ckq(i)  can  he 
written  in  terms  of  those  tensors  as 


and 


where 


Ckq(i)  =  /2^M  <£!ICklli>u*(i) 


l  Ckq(i)  «  m+1  <UCkll?,>U^k) 


'k  q 


(115) 


(116) 


<JUC  lli>  =  <£(0)k(0)  |  l{  0)>  . 

K 


Thus  equation  (114)  may  be  written 


H3  =  k|  <*«Ck!U>U<k>  . 


(117) 
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The  reduced  matrix  elements  tabulated  by  Neilson  and  Koster 
are  for  L-S  states,  so  we  must  relate  the  matrix  elements  in  the  J  basis 
to  L-S  basis.  To  do  this  we  construct  the  state 


Then 


|  JMLS>  *  V  <L(p)S(M-p) |J(M)> |LpSM-p> 
P 


(1*8) 


CJ'M'L'S'  (l/k)  |JMLS>  ■  <J(M)k(q) | J’ (M* )><J'L'S' DU*k) KJLS>  ,  ,  (119) 

by  application  of  the  Wigner-Echart  theorem.  If  now  we  calculate  the 
matrix  element  in  equation  (119)  by  using  the  wavefunctions  of  Equation 
(118),  we  have 

<J'M,L'S' |U(k) | JMLS> 

q 

=  l  <L(  p)S(M-p)  |J(M)><L'  (  jj  *  )S'  (M*-p'  )  |J'  (M*  )> 

UP’  (120) 

x  <L' p'S'M'-p' |U<k) |LpSM-p>  . 

q 


The  matrix  elements, 

( k  1 

<l' u's'm'-m1 |ir  '|LpSM-p> 

*  6SS'6M'-u',M-y<L(y)k(q),L’U!,)><L,S,,lu(k)"LS>  ' 


(  121) 


are  obtained  using  the  Wigner-Echart  theorem,  with  the  knowledge  that 
is  a  spherical  tensor  in  orbital  space  only.  In  equation  (120)  the 
last  C-G  coefficient  and  the  C-G  in  equation  (121)  can  be  rearranged  to 
give 


<L(  p)k(q)  |L‘  (p'lXV  (  U'  )S(M-p)  |  J'(M,)> 

T  4-lr—T  * 

=  (-1)  <k(q)L(p)|L'(p')XL'(p*)S(M-p)|j'(M,)>  , 


where  we  have  used  the  results  S'  «  S  and  M'  -  jj  *  *  M  -  p  given  in 
equation  (121).  The  two  C-G  coefficients  given  in  equation  (122)  can  be 
written 
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L-fk-L  * 

(-1)  <k(q)L(p)  |LMU' }><L*  (U*  )  S(M-|i )  |  J '  ( M'  )  > 

-  (-1)L+k"L  l  /?2L'+1)(2f+1)  W(kLJ'S»L'f )<L(m)S(M-M) |f(M)>  (123) 

f 

x  <k(q)f (M) | J' (M* )>  . 


Since  all  the  dependance  on  y  is  contained  in  the  first  C-G  in  equation 
(120)  (the  only  one  remaining)  and  the  first  of  equation  (123),  this  sum 
can  be  performed  to  give 

l  <L(w)S(M-y)  |j(M)XL(y)S(M-y)  if  (y)>  *  6Jf  (124) 

U 

from  the  orthogonality  of  the  C-G  coefficients.  Substituting  these 
results  back  into  equation  (120)  we  have 

CJ'M'L’S'  |U^k)  |JMLS>  =  (  - 1 ) L+k-1< '  /(  2L '  + 1  )  (  2  J+ 1 )  W(kbJ' Sf  L' J) 

x  <k(q)J(M)  i  J*  (y1  )><L'SIIU(k)  HLS>6SS,  .  (125) 

We  have  used  the  condition  f  =  J  in  equation  (124)  and  the  result  in 
equation  (121).  To  compare  the  result  in  equation  (125)  with  equation 
(118),  we  need  to  rearrange  the  C-G  in  equation  (125): 


<J'L'Sl1U<k)  BJLS>  =  (  -1  )L+J'  ~'L,~  J/('2L,  +  1 )  (  2J+1 ")  W(kLJ'S;L'J) 

x  <L'SIU(k) «LS>  . 


(  126) 


The  matrix  elements  tabulated  by  Neilson  and  Koster  (1964)  are  related 
to  ours  by 


(li'SIU(k)  HLs)  =  /Sl'  +  I  <L'SIIU(k)  BLS> 


(  127) 


Finally  the  complete  matrix  elements  of  H3  given  in  equation  (117)  are 
<j'm'l's|h3|jmls> 

£  BiL<J(M)k(q) |J’ (M' )>  /2l+T  <AICtli> 
kq  Kq  k  (128) 

x  <J'L'SHU(k)IJLS>  , 
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with 

<J'L’SI!U(k)  itJLS>  "  (-1)L“L'+J,-J/2JTT  W(kLJ'S;L'J)  (129) 

x  (L'SHU(k)  II LS  ) 

and 


<£IICk#£>  =  <«{0)k(0)  U(0)>  . 

For  i  =  3,  < £ HCk « «,>  =  -2//T5 ,  2//22,  and  -10//429  for  k  =  2,  4,  and  6, 
respectively. 

4  3+ 

4.1.2  Numerical  Example:  Fj  States  of  Nd 


As  a  numerical  example  of  the  calculation  of  the  crystal- 
field  splitting,  we  will  calculate  the  splitting  for  Nd3+  (4f3)  in  a 
field  of  S4  point  symmetry.  We  shall  assume  that  the  levels  are  pure 
^Fj.  We  assume  that  L,  S,  and  J  are  all  good  quantum  numbers;  then  we 
consider  matrix  elements  of  in  equation  (128)  with  J'  =  J  and  L'  = 
L.  Thus, 


<JM'LS|H3I JMLS>  =  l  B*  <J(M)k(q) | (M' )> 

kq  (130) 

x  <JIIUkIIJ>  /2l+T  <£(0)k(0)  U(0)>  . 


The  values  of  the  reduced  matrix  elaments  <JHUkIIJ>  /2l+1< It  110^11  £>  for  the 
V  state  of  Nd3+  are  as  in  table  9.  In  obtaining  these  values  we  have 
used  Neilson  and  Koster's  results  (1964)  for  the  reduced  matrix  ele¬ 
ments,  (LSllUkliLS) ,  for  L  =  3  and  S  —  3/2;  the  Racah  coefficients  (6 
symbols)  are  found  in  Rotenberg  et  al  (1969). 

TABLE  9.  VALUES  OF  <  J  >Uk  I  J>  /2l+T<  i  DC^  II  l> 


J  k  =  2  k  =■  4  k  =  6 


3 

2 

5 

2 

7 

2 

9 

2 


\  & 


1/2 


UL-  r_i  1 

2*5  12  *5  J 


1/2 


3  [3^] 


To  ijo] 


j/2 

1/2 


1  r__L_ 

12 -3  • 


r 

,1/2 


2  L  2  •  3  « 

_l_  r_ia_|l/z 

2.3  1-2  - 1 1  j 


0 

0 


,l/2 


_  I2I  r _ 3 _ 1 

3  I3  »1 1  *13  J 

i  r _ 1 _ 1^2 

3  L2-3-11 *13J 


The  calculation  of  the  energy  levels  is  made  somewhat  simpler 
by  using  wave  functions  that  transform  according  to  some  irreducible 
representation  of  the  S4  group  (Koater  et  al,  1963).  The  irreducible 
representations  of  the  S4  group  are  all  one-dimensional,  but,  since  the 
ion  under  investigation  has  an  odd  number  of  electrons,  the  energy 
levels  will  be  at  least  doubly  degenerate.  Thus,  of  the  four  irreduc¬ 
ible  representations — Fg,  Fg,  Tj,  and  Tg — only  two  need  be  chosen;  we 
chose  Fg  and  F-^.  The  energy  corresponding  to  Tg  is  degenerate  with  Tg 
and  that  corresponding  to  rQ  is  degenerate  with  The  wave  functions 

belonging  to  T-,  with  a  particular  J  value  are 

U  a|±j  >  ,  _  mi  <  k  < 

and  those  for  are 

u  ^  ,  _  an  <  k  <  ma  , 

where  k  is  an  integer,  and  the  number  of  k  values  occurring  for  a  given 
J  is  the  number  of  times  a  representation  will  occur.  The  number  of  3kq 
for  the  calcium  site  in  calcium  tungstate  is  five:  i.e.,  b20,  B4Q,  B60' 
B44,  and  Bg4 •  Of  these  parameters  only  the  Bg4  is  complex.  The  matrix 
elements  of  the  crystal  field  given  in  the  above  equation  are  presented 
explicitly  under  particular  states  in  the  following  paragraphs. 


This  level  of  the  free  ion  is  split  into  two  doublets  by  the 

crystalline  field.  The  wave  functions  corresponding  to  and  to  T,. 
_  I  3  1  V _ J  I  3  3  k  _ _  1  1 _ 1  mm _  _ . _ /  *  \  I _  . 


are  j  >  and  1^  ~  'J  >  respectively.  From  equation  (130)  we  have 

(  H  ^  H3 1  W  >  =  ”  2?  B20  =  E^r7  2^ 


(  131) 


<!-!iHj!-!>=fgB2o  =  E(F,  |) 


(132) 


where  the  appropriate  values  of  the  reduced  matrix  elements  in  equations 
(131)  and  (132)  were  taken  from  table  9.  The  total  splitting  of  the 
**F3/2  state  *-s  then 


—  b 
25  20 


(133) 
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5/2 


Unlike  the  previous  case,  this  state  contains  two  T.'s,  and 
their  wave  functions  are  1^  2  >  and  1^  “  ^  >•  The  wave  function  for 
the  r7  state  is  I'f  ^  The  energy  for  r7  is 


<  H  |H3 1  f  5  >  -  7(10  f"4B20  *  If  B4qJ  ■  E(r7  f)  •  (’341 


5  1 


11 
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The  necessary  matrix  elements  for  the  energy  in  are 
.  5  5  ...  ,  5  5  .  11 


<5.3 
2  2 


5  3 


t5B20 

4  38  B40^  b11  ' 

(135) 

-11  r, 

700  L 

"B20  "  TT  B40  ^  =  b22  ' 

(  136) 

<f-f  '"3'H>--=»44-b,2 


1/2 


The  two  energy  levels  corresponding  to  r5  are 

F  fr  S'!  _  bn  +  b22  *  i(b11  ~  b22 ) 2  +  4b  12bt 2 ] 

EUI5  2J  ~  2 

v  (r  1)  _  b11  4  b22  ~  [(b11  ~  b22 ) 2  *  4b12b^2]l/2 

E2^5  2>  ~  2 


(137) 


(138) 

(139) 


7/2 


This  state  contains  two  T^'s  and  two  Ty's  . 

elements  for  Tj  are 

v  7  1  .  .  7  1  1  r  55  n  ,  9  _  ,  50  _  1 

<  2  2  2  2  >  =  99  *■”  7  B20  14  ®40  13  B60  J  *11  ' 

<  2  ”  2  'H3'  2~2>  99  [11B20  +  2  B40  "l3  B60  ^  =  ®22  ' 


<1-1  lH3l  \  \> 


_2 

99 


-1b  -  2£  /io  b 

^  44  13  64 


‘12 


The  matrix 

( 140) 

(141) 

(142) 
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The  two  energy  levels  are 

E1(r?  J)  -  a11  +  a22  *  Ka11 
E2(r7  |)  -  *)J  *a-?2  ~-[(a” 

The  matrix  elements  for  are* 


~  a22)2  *  4a 12a^2 ] 1/2 
2  ' 

"  a22)2  +  4a12a!2JM/2 

T. - * 


(143) 

(144) 


<  1  f  |H3 

1 

1  2  2 

>  «-L 

99 

B20 

-  22  B 

14  40 

+  sa 

13 

B60^  ”  b11  ' 

(  145) 

<  l 

|H3'  \ 

_  2  > 

1 

I"  ”  *20 

- —  B 

14  40 

90 

B6q]  =  b22  ' 

(146) 

<  2 

"  99 

13 

< 

\~\  |H3' 

21, 
2  2 

1 

99 

Vri  ®44 

^  20  I21 
+  1*3  \~2 

B64] 

=  b12  * 

(147) 

The  corresponding  energies  are  given  by  substituting  the  above  values  of 
b^j  into  equations  (138)  and  (139). 


The  number  of  r7's  in  this  state  is  three,  with 
The  matrix  elements  of  the  crystal  field  for  r7  are 

<  ?  ?  |h3!  \  f  >  =  [6B20  "  ^  ®40  +  ^  B6o]  =  a11  ' 

<f^lH3lf^>“  3^  [_4B20  "  ^  B40  -  B6oJ  1  a22  • 

<  I  "  J  |H3!  \  "  \  >  "  3^5  [2B20  +  B40  "  ^  B6o]  "  a33  ' 

<  f  i  lH3l  f  I  >  "  ife  [■  f  ^  B44  +  ^  B64  ]  "  a12  ' 


two  r5’s. 

(148) 

(149) 

(150) 

(151) 


<  f  "  i  |H3'  I  i  >  “  ‘sfe  ["  A  B44  +  B64]  ■  a23  *  (152) 


*The  symbol  a±.  will  be  used  for  the  matrix  elements  in  ?7  and  b±j 
for  those  in  r5  to  avoid  introducing  new  symbols  for  each  new  value  of 
J. 
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The  three  energies  are  given  by  the  solutions  of 


E3  -  (a1 1+a22+a33)E2  +  (a  i 1a22+a 1 1a33+a22a33-a23a23~a 12M2 )E 
+  a11a23a:$3  +  a33a 1 2a1 2  “  a11a22a33  =  0 


The  matrix  elements  for  Tg  are 


(  153) 


2  2  3  2  2 

2  _  3  |H  I  2  _  3 

2  2  3  2  2 


>  "356 

[B  +  12  B  +Mb 

L  20  7  40  91  60 

(  154) 

>  "Be 

[”3B20  “  7  B40  +  97  B60 

^  =  b22  ' 

(  155) 

>  .  J L 

396 

t-  f  /5°  S44  '  IT  /6  B64 )  *  b,2  • 

(  156) 

and  E2(r,j  2 )  are  9^ven  by 
using  the  b^j  given  above. 

equations 

(138)  and 

Calculation 


I 

t 

i 

i 

j 

I 

i 


Three  of  the  crystal-field  parameters  can  be  obtained  quite 
simply  from  the  experimental  data.  These  are  B2g,  B4Q,  and  B^g.  if  we 
express  the  sums,  S^,  in  terms  of  the  where  the  are  experimental 
data,  then  we  obtain 


and 


SJ-1/2  =  l  Ei(r7J) 


S1  "  “  25  ®20  ' 


(157) 


„  11  ,  1  „ 

S2  175  ®20  42  B40 


158) 


_2  40 

63  20  693  40  1287  60 


(  159) 


_  _Z_  B  _  _40_ 

4  9  20  198  40  1287  60 


(160) 


where  S4  is  for  the  J  =*  9/2  level.  These  equations  can  be  inverted  to 
give 


B 


20 


-  21  S 

2  S1  ' 


(161) 
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:  1 

i  i 


B40  =  "33S1  +  42S2 


(162) 


B  _  JL20J.  _  78  J287 

B60  40  S1  5  S2  +  40  S3 


(163) 


The  other  crystal-field  parameters  are  slightly  more  in¬ 
volved.  Prom  equations  (138)  and  (139)  we  have 


where 


=  _  11 


B44  -  3^4  [wf  -  Nf]l/2  , 


S,  +  2S2  and  W,  =  f)  -  E2(r5  f)  . 


(164) 


To  determine  Bg4,  we  use  equations  (143)  and  (144)  to  give 

a12a12  =  4  lW2  ~  N2 ] 


(165) 


where 


"2  ■  5  S1  -  f  S2  *  I  S3  *"d  «2  ■  El(r7  1)  -  B2(r7  1)  • 

A  similar  expression  can  be  obtained  using  equations  (145) 
and  (146)  in  equations  (138)  and  (139),  yielding 


b12b12  -  4  tW3  “  ^3 ]  ' 


(166) 


where 


“3  '  S1  '  2S2  *  1  S3  *"d  W3  '  El(F5  ?)  -  E2^F5  • 

Substituting  equations  (142)  and  (143)  into  the  left  side  of 
equations  (165)  and  (166),  we  obtain  two  equations  for  Bg4.  These  two 
equations  can  be  solved  simultaneously  for  both  real  and  imaginary  parts 
of  Bg4  to  give 

*6*4$  ft22  (W2  +  W3  “  KI  “  N§)  -  6(W1  -  Nl)]l/2  '  (167> 


cos  0  = 


1 5  *99 
13  8 


(W§  -  M|)  -  8(w2  -  N2)  -  1^-  (W2  -  N2) 


20/70 


Vi  -  »?)1/2 


1168) 


where 
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All  the  crystal-field  parameters  can  be  determined  once  the 
experimental  data  are  taken  on  the  '+F5/2'  an<J  4p7/2  ]Lavell5• 

As  tedious  as  the  above  procedures  may  have  seemed,  the 
crystal -field  parameters  we  obtain  are  only  approximate  since  we  have 
ignored  L-S  mixing  by  the  spin-orbit  coupling  in  the  free  ion  and  J 
mixing  caused  by  the  crystal  field.  Nevertheless,  the  crystal-field 
parameters  obtained  by  the  above  procedure  can  serve  as  very  good  start¬ 
ing  values  in  a  fitting  of  a  more  sophisticated  calculation  to  experi¬ 
mental  crystal-field  levels. 

The  crystal-field  parameters  B^g  obtained  by  the  above  proce¬ 
dure  for  Nd3+  in  CaWO^  are  given  below,  along  with  crystal-field  param¬ 
eters  for  the  same  ion  but  with  full  diagonalization,  that  is,  L-S 
mixing  and  J  mixing  (Wortman  et  al,  1977). 


B20 

B40 

®44 

B60 

RB64 

IB64 

above 

403 

-635 

+711 

-219 

885 

0 

full  diagon- 

509 

-866 

1042 

-509 

903 

243 

alization 

4 . 2  Classical  Point-Charge  Model 


In  the  simplest  model  of  the  crystal  field,  the  point-charge 
model  introduced  by  Bethe  (1929),  the  lattice  is  replaced  by  an  array  of 
point  charges  placed  at  the  nuclei  of  the  constituent  ions.  A  multipole 
expansion  is  made  of  the  point-charge  potential  energy  at  the  rare-earth 
site.  Thus,  if  R^^ ( j)  is  the  vector  position  of  constituent  j  at  site 
j  in  the  i,  m,  nth  cell,  we  have 


H 


3 


l  l 

imn  j 


(j)  -  r| 


(  169) 


where 


+  ♦ 

Jta  +  mb  +  nc  + 


+• 


>  >  >  -*■ 
and  a,  b,  and  c  are  lattice  vectors.  The  charge  at  site  j  is  eZj  and  r 

is  the  position  of  an  electron  on  a  rare-earth  ion.  The  multipolar 

expansion  of  equation  (169)  is 
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(  170) 


-e2z,rk 


C  j  i-  ^ 

'fin  j  [R£m<  3)  ]k+1  Ckq<r)Ckq[^mn(^]  * 


The  multipolar  crystal-field  components  A^g  are 


Akq  = 


I  I 


tan 


2jCkqLiW3>] 

[Rimn^>]k+1 


(171) 


Thus  the  point-charge  Hamiltonian  is 


N 

H3=IAkq  l  riCkq(^i)  '  <172> 

kq  i  =  1 

N 

where  we  have  summed  over  all  the  N  electrons  in  the  4f  configuration. 
If  all  the  lengths  are  measured  in  A  (10~8  cm),  then 


x  10 


8 


1 

Amn 


I  r 

j  [RAmn(j)] 


k+1 


( 173) 


where  aQ  =  e2/tfc  (a0/2ir  x  108 
constant — and  the  units  of  A^g 
Akq<rk>  =  cm'1. 


1 1  6, 1 40)-~that  is,  the  fine-structure 
cm'1/^*  If  <rk>  is  in  A  units,  then 


The  sum  in  equation  (173)  always  converges — even  for  the  lowest 
k  value  (k  *  0)— if  taken  in  the  order  indicated.  That  is,  the  sum  over 
j  is  performed  with  A,  m,  and  n  fixed.  The  unit  cell  is  neutral,  that 
is. 


)  Z.  =  0  .  (174) 

j  j 

In  many  cases  (not  all  space  groups)  it  is  possible  to  choose  an  origin 
for  the  lattice  coordinates  such  that  the  dipole  moment  of  the  unit  cell 
vanishes;  that  is. 


Pjzj  -  0  ,  (175) 

where  is  the  position  of  the  jth  charge  in  the  unit  cell.  The  result 
in  equation  (175)  can  be  anticipated  by  observing  the  point  symmetry  of 
tho  ions  in  a  specific  solid.  If  the  ions  occupy  C^,  C2»  Cg,  C2v,  C 4, 
C4v'  C3'  C3v'  C6'  or  C6v  Point  symmetry  (Schoenflies  notation),  then  it 
is  impossible  to  satisfy  equation  (175)  with  these  sites  as  the  origin 
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in  a  unit  cell.  If  it  is  possible  to  satisfy  equation  (175).  then  the 
sum  given  in  equation  (173)  converges  very  rapidly.  This  can  be  shown 
from  the  expansion  (Carlson  and  Rushbrooke,  1950) 


l«-l!kTl 


<a(  a)k(q)  |a+k(  oH-q)>x3C  (x) 

ci  ex 


c* 

a+k , g+g 
Ra+k+1 


/s 

(R) 


(  176) 


With  =  x  and  $£mn(0)  =  $  n(j)  =  ^jjmn(0)  +  Pj)k  for  the  sum  in 

equation  (173)  we  have 


7  7  Z jCkq ^ R£mn^  ^  ^  7  7  />2a+2k\1'/2 

i  l  ~T  a - 7J—  =  l  1  V  2a  j  <a  (  a)k  (q)  |  a+k  (  n+q )  > 

£mn  aa  '  ( 1 


77) 


> !  !.*  (;.) 
j  3  3  *“  3  [r  (0,l*tk+1 

L  £mn 


(0)  ] 


Now  if  equation  (175)  is  satisfied,  then 

T  Z.p.C,  fp.)  =  0  .  ( 178) 

j  3  ^  1a'  3 

Thus,  the  sum  in  equation  (177)  is  for  a  >  1 ?  we  see  that  even  for  the 
lowest  term,  k  =  0,  we  have  the  individual  terms  falling  off  as 
1/R|mn(0).  While  the  expansions  in  equations  (176)  and  (177)  are  good 
for  demonstrating  the  rate  of  convergence,  the  computation  of  A,^  by 
equation  (173)  is  more  practical.  However,  in  equation  (173),  the  sum 
over  j  should  be  done  for  each  cell  first,  with  fixed  values  of  £,  m, 
and  n.  In  programming  language,  this  is  expressed  by  stating  that  the  j 
loop  is  the  innermost  of  the  nested  £,  m,  n,  and  j  loops.  In  some 
lattices,  the  condition  in  equation  (175)  may  place  some  of  the  point 
charges  on  the  cell  faces.  In  these  cases  it  is  a  simple  matter  to 
balance  these  charges  by  an  adjustment  to  fractions  of  equal  charges  on 
opposite  faces. 


The  convention  we  use  for  our  lattice  sums  is  that  given  in  the 
International  Tables  for  Crystallography  (1952);  table  10  is  reproduced 
from  volume  I  (the  other  two  volumes  give  data  strictly  for  x-ray  crys- 
tallographers ) .  The  data  used  in  the  lattice  sums  are  generally  those 
reported  in  Acta  Crystallographies,  Section  B.  Care  should  be  taken  to 
make  certain  the  correct  setting  is  used. 


Typical  data  used  in  the  calculation  of  the  lattice  stuns  are 
given  in  table  11  (LiYF^,  calcium  tungstate  space  group  88,  Scheelite 
structure).  All  the  data  given  in  table  11  are  given  in  the  Interna¬ 
tional  Tables,  except  that  the  x,  y,  and  z  coordinates  are  determined  by 
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TABLE  10.  CRYSTALLOGRAPHIC  AXIAL  AND  ANGULAR  RELATIONSHIPS  AND 
CHARACTERISTIC  SYMMETRY  OF  CRYSTAL  SYSTEMS 


Space 

group 

System 

Axial  and  angular 
relationships 

X-ray  data  needed 
for  unit  cell 

1,2 

Triclinic 

a  *  b  *  c 
a  *  0  *  y  t  90° 

a,  b,  c,  o,  B,  Y 

3  to  15 

Monoclinic 

First  setting 
a  *  b  *  c 
a  =>  6  =  90°  *  y 

Second  setting 
a  *  b  *  c 
a  =  y  =  90°  t  6 

a,  b,  c,  y 

a,  b,  c,  B 

16  to  74 

Orthorhombic 

a  *  b  *  c 
<x  =  3  =  y  =  90° 

a,  b,  c 

75  to  142 

Tetragonal 

a  =  b  *  c 
a  =  B  =  Y  =  90° 

a ,  c 

143  to  167 

Trigonal 
(may  be  taken 
as  subdivision 
hexagonal ) 

(Rhombohedral  axes) 
a  =  b  =  c 

a  =  6  =  Y  <  120°  *  90° 
a  =  b  *  c 
a  =  6  =  90° 

Y  =  120° 

a ,  a 

168  to  194 

Hexagonal 

a  =  b  *  c 
a  =  6  =  90° 

Y  =  120° 

A,  C 

195  to  230 

Cubic 

a  =  b  =  c 
a  =  B  =  Y  =  90° 

a 

Source:  International  Tables,  1952,  Vol ,  I,  p.  11,  table  2,3.1. 


x-ray  diffraction.  The  lattice  constants  a,  b,  and  c  are  also  deter¬ 
mined  by  x-ray  diffraction  and,  as  customary,  the  true  positions  of  the 
ions  are  xa,  yb,  and  zc  (these  relations  hold  for  all  the  ions  in  a  unit 
cell).  The  polarizability  (from  Kittel,  1956;  165)  of  each  ion  is  given 
at  the  bottom  of  table  11.  For  this  particular  solid,  only  the  fluorine 
ions  can  have  dipole  moments  that  contribute  to  the  crystal  field  (we 
shall  discuss  the  role  of  the  dipole  moments  later  on).  All  the  data 
for  space  group  88  are  not  contained  in  table  11  because  the  equivalent 
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positions  given  in  the  International  Tables  are  generated  inside  the 
program.  The  centering  position  in  the  cell  can  be  taken  as  either  the 
Y  or  Li  site,  since  these  positions  have  S4  symmetry,  and  their  lowest 
crystal-field  component  is  Ajg.  Equation  (178)  is  therefore  automati¬ 
cally  satisfied.  The  resulting  lattice  sum  for  the  Y  site  in  LiYF^  for 
the  parameters  in  table  1 1  is  given  in  table  12. 

TABLE  11.  CRYSTALLOGRAPHIC  DATA  FOR  LiYF4  (SCHEELITE,  CaW04) 


TETRAGONAL 

SPACE  GROUP 

08  (FIRST 

SETTING)  Z 

=*  4 

Ion 

Position 

Symmetry 

X 

y 

z 

Y 

4b 

S4 

0 

0 

1/2 

Li 

4a 

S4 

0 

0 

0 

F 

16f 

C1 

0.2820 

0.1642 

0.0815 

Notes:  a  =•  5.1668,  b  -  a ,  c  =  10.733,  a  *  9C°,  0  =  90°, 
Y  =  90°/  a  »  0.55  A3,  aLi  ■  0.05  A3,  aF  =  1.04  A3  (reduced 
to  0.104  in  the  lattice  sum),  Zy  =  +3,  ZLi  =  +1,  Zp  =  -1. 


TABLE  12.  LATTICE  SUMS  FOR  Y  SITE  AT  (0,  0,  1/2)  FOR 
LiYF4  WITH  Zy  *  3,  -  1  ,  Zp  ■  -1 ,  *  0.104  A3 


Lattice 

sum 

Monopole  Anm 

Dipole  Anm 

Monopole  and  dipole 

Real 

Imaginary 

Real 

Imaginary 

Real 

Imaginary 

A20 

1074 

0 

340 

0 

1  414 

0 

A32 

373 

859 

-358 

74.0 

15 

933 

a40 

-1957 

0 

-98.1 

0 

-2055 

0 

a44 

-2469 

-2362 

-3.83 

-80.3 

-2473 

-2442 

A52 

1050 

-2456 

1.28 

-74.7 

1051 

-2531 

A60 

-17. 

2  0 

7.96 

0 

-9.24 

0 

a64 

-615 

-420 

-29.03 

-9.  37 

-644 

-429 

A72 

-15. 

7  0.90 

1.55 

-9.94 

-14.2 

-9.04 

A76 

250 

-63.8 

1  7.8 

7.96 

268 

-55.9 

The  sum  covers  all  of  the  complete  cells  in  a  sphere  of  30- 
A  radius  and  should  be  an  accurate  result.  Also  included  in  table  12 
are  the  results  for  the  dipole  contributions  due  to  the  presence  of 
dipoles  at  the  fluorine  sites  (we  shall  discuss  these  terms  later). 
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ppm  jMin  pun  ill 


As  a  second  example,  we  choose  a  very  low  symmetry  solid,  YCl-j, 
which  is  monoclinic  space  group  12.  As  can  be  seen  in  table  13,  all  the 
ions  are  in  very  low  point-symmetry  positions,  and  each  position  can 
have  a  dipole  moment  (another  way  of  saying  this  is  that  each  position 
has  a  onefold  field,  A1m) .  We  then  have  to  consult  the  International 
Tables  for  a  higher  symmetry  position  in  order  to  satisfy  equation 
(175),  which  in  this  case  is  the  site  4e  with  symmetry.  The  point 
group  has  only  the  inversion  operation,  and  all  the  odd-n  An  vanish  in 
this  symmetry.  Thus  if  the  position  4e  is  used,  equation  (175)  will 
automatically  be  satisfied.  The  lattice  sum  for  YC13  was  also  run  over 
a  lattice  30  x  30  x  30  and  only  the  even-n  Anm  are  given  in  table  14. 
The  dipole  contributions  were  also  calculated;  these  calculations  were 
more  complicated  in  this  solid  because  of  the  three  types  of  sites  (Y, 
Cl^,  Cl2).  For  many  of  the  A^,  the  dipole  contributions  are  much 
larger  than  the  monopole  terms.  This  frequently  happens  when  the  hand¬ 
book  values  for  the  dipole  polarizabilities  are  used.  We  have  had  more 
believable  results  when  we  reduce  the  polarizability  to  one  tenth  of  the 
handbook  value. 


TABLE  13.  CRYSTALLOGRAPHIC  DATA  FOR  YClj 
MONOCLINIC  SPACE  GROUP  12  (C2/m)  (SECOND 
SETTING)  Z  -  4 


Ion 

Position 

Symmetry 

X 

y 

z 

Y 

4g 

C2 

0 

0.166 

0 

cii 

Cl2 

4i 

cs 

0.211 

0 

0.247 

83 

C1 

0.229 

0.179 

0.760 

4e 

ci 

1/4 

1/4 

0 

Notes : 

<3  *  S  • 

92,  b  *■  11 

.94,  c  -  6. 

44,  a  “ 

90,  (3  - 

111.0, 

y  -  90. 

Charges:  gy  »  3,  qc^  »  -1. 

Polarizability:  &y  «  0.55  A3,  *  3.66  A3. 


TABLE  14.  LATTICE  SUMS  FOR  Y  SITE  AT  (0,  0.166,  0) 
FOR  YClj,  EVEN-n  Anm  ONLY,  ALL  Anm  REAL 


Lattice  sum 

Monopole 

Dipole 

Total 

A20 

1738 

3227 

4965 

a2  1 

-913 

2916 

2003 

A22 

245 

2574 

2819 

a40 

-73.9 

246 

172 

A4 1 

85.8 

-398 

-312 

a42 

-41.3 

47.7 

6.4 

A4  3 

10.4 

-791 

-781 

a44 

-3.64 

516 

512 

A60 

-0.06 

-80.2 

cn 

o 

CD 

1 

A61 

-3.76 

21.3 

17.5 

A62 

3.35 

-2  /.4 

-24.0 

A6  3 

-0.58 

60.5 

59.9 

A64 

2.73 

39.2 

41.9 

A65 

5.07 

13.7 

18.8 

A66 

8.  14 

-65.9 

-57.8 

The  lattice  sums  given  in  tables  12  and  14  are  incomplete  in 
that  the  results  are  not  in  a  usable  form  for  many  of  our  computer 
programs.  Before  we  can  use  these  results,  the  Anm  should  be  rotated 
using  the  standard  Euler  angle-rotation  matrix,  so  that  the  lattice 
sums,  A^,  rotated  from  Anm  by  the  angles  a,  3,  and  y  are 


A' 

nm 


0  (a,3,y)A 

m  m  nm 


Explicit  forms  for  the  Dn,  (a,  0,y)  are  given  in  Rose  (1957,  ch.  IV). 

mm 

4.3  Point-Charge  Model  Developed  at  HDL 
4.3.1  Introduction 


In  the  classical  point-charge  model,  the  crystal-field  param¬ 
eters,  Bnm,  for  the  crystal-field  interaction  of  the  form 


were  calculated  as 


-X 


B*  l 

nm  ~ 
nm  i 


cm 

nm 


B  =  <rn>A  , 
nm  nm 


(  179) 


(180) 


where  the  <rn>  are  the  expectation  values  of  rn  of  the  rare-earth  ion, 
and  the  Anm  are  the  multipole  components  of  the  energy  at  the  site 
occupied  by  the  rare-earth  ion.  In  the  earlier  models,  the  radial 
integrals  used  in  the  evaluation  of  rn  were  taken  from  Hartree-Btock 
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calculations  (Freeman  and  Watson,  1962),  and  the  were  calculated 
using  the  point  charges  at  the  valence  values  tor  the  constituent  ions. 
These  calculations  generally  gave  the  twofold  field  10  timus  too  large, 
the  fourfold  fields  approximately  in  good  agreement,  and  the  sixfold 
fields  10  times  too  small. 

4.3.2  Screening  and  Wave  Function  Spread 

Several  errors  in  the  classical  theory  were  immediately 
obvious.  If  the  radial  wave  functions  (Hartree-Fock)  for  the  free  ion 
were  correct,  then  these  wave  functions  should  give  the  correct  values 
for  the  Slater  integrals  F**,  F4,  and  F®.  They  did  not  for  Pr^  +  .  A 
simple  procedure  was  then  applied.  The  radial  wave  functions  were 
assumed  to  be  of  the  form 


♦  <r)  =  CRHp(Tr) 


(  181) 


where  r  is  a  parameter,  C  is  a  normalization  factor,  and  R^fr)  are  the 
Hartree-Fock  radial  wave  functions.  With  the  radial  function  given  by 
equation  (181),  the  Slater  integrals  become 


F 


k 


(  182) 


and  it  was  found  that  a  value  of  r  of  approximately  0.75  was  needed  to 
fit  the  Fk  that  are  found  by  fitting  the  experimental  spectra  of  Pr3f. 
Thus,  the  Hartree-Fock  radial  wave  functions  had  their  maxima  too  near 
the  origin  and  needed  to  be  spread  out  even  in  the  free  ion,  and  perhaps 
more  spreading  was  necessary  for  an  ion  in  a  solid. 


From  the  radial  wave  functions  given  in  equation  (181),  it  is 
not  difficult  to  show  that 


<f(r)>  *  <4>  |  f  (  r )  l4>>/<4>|<f>> 


=  <f(r/T)> 


HF 


(  183) 


so  that  any  quantity  that  has  been  calculated  using  Hartree-Fock  func¬ 
tions  is  immediately  obtained,  particularly 

<rk>  =  <rk>HF/xk  .  (184) 

A  second  error  of  the  classical  method  was  the  omission  of 
the  Sternheimer  shielding  factors  ( Sternheimer,  1951,  1966?  Sternheimer 
et  al,  1968).  In  1951  Sternheimer  showed  that,  in  a  multipolar  expan¬ 
sion  of  the  energy  of  a  point  charge  embedded  in  a  solid,  the  rn  should 
be  replaced  by  rn(l  -  onJ,  where  the  an  are  the  shielding  factors.  He 
further  showed  that  these  factors  were11  independent  of  azimuthal  angle? 
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that  is,  if  the  angular  variation  in  the  multipolar  expansion  was  given 
by  Y  ,  the  on  were  independent  of  m.  The  values  of  0  have  beer,  calcu¬ 
lated  for  Pr3+  and  Tm8+  and  are  1 

3+ 

a.  *  0.666,  a.  *  0.09,  a,  -  0.04  for  Pr  , 

2  4  b 

3-  (185) 

a_  «  0.545,  a.  »  0.09,  a.  m  0.04  for  Tm  , 

2  4  6 

where  the  replacement  is 

rn  -*•  rll(  1  -  a  )  .  ( 186) 

More  recent  calculations  of  the  shielding  factors  have  been  done 
(Sengupta  and  Artman,  1970,  and  perhaps  others),  which  we  shall  need  if 
further  refinements  of  the  theory  are  undertaken. 


4.3.3  Effective  Charge  and  Position 


The  crystal-field  components,  Anm,  are  a  function  of  the 
position  of  the  ions  in  a  solid;  in  solids  such  as  CaWO^  tne  (WC/^)”2 
complex  is  known  to  be  covalent.  That  is,  the  charges  on  the  tungsten 
and  the  oxygen  ions  are  not  necessarily  at  their  valence  values.  If  we 
let  the  charge  on  the  tungsten  ion  be  q^,  then  we  require  that 

qW  +  4qo  ”  "2  (187) 

with  the  charge  on  the  oxygen  being  q^.  The  result  given  in  equation 
(187)  then  assumes  that  the  Ca2+  site  is  purely  ionic  with  charge  2. 
This  assumption  is  consistent  with  many  of  the  experimental  results  on 
compounds  such  as  CaWO^  or  YVC>4.  We  introduced  a  second  parameter,  the 
effective  position  of  the  oxygen  ion  relative  to  the  tungsten  site  that 
would  reproduce  the  effective  dipole  moment  seen  from  the  Ca  site. 
This  parameter,  rt,  is  introduced  by 


R  (effective)  =  nR,,  (measured)  ,  (188) 

U*W  w“W 

where  Rq_w  is  the  distance  from  the  oxygen  nucleus  to  the  tungsten 
nucleus.  Thus  there  are  two  parameters  in  the  ^  ;  q^,  the  effective 
charge  (1^  is  eliminated  by  eq  (1873),  and  n»  the  effective  distance  of 
the  oxygen  site  from  the  tungsten  site.  The  calculated  crystal-field 
parameters  Bnm  then  are 


Bnm(T;a011)  "  <rn>Hl?(1  "  ojArun^O'^” 


with  the  three  parameters  t,  qQ,  and 


( 189) 


56 


"n^w-nri 


The  experimental  data  that  were  taken  at  HDL  on  the  rare-earth 
ions  in  CaWO^  were  analyzed  using  the  effective  spin-orbit  Hamiltonian 
(Karayianis,  1970),  and  a  set  of  phenomenological  Bnm  was  obtained. 
These,  given  in  table  IS,  are  the  Bnm  that  the  theory  has  to  fit. 

The  fitting  was  done  by  minimizing  the  square  quantity  given  by 

Q  “  HBnm('‘»qon)  -  BmJ2  '  <190> 

nm 

where  9iven  by  equation  (189),  and  Bnm  is  from  table  15, 

for  each  ion.  The  minimization  was  done  with  respect  to  r,  qQ,  and  n 
for  each  ion.  Since  the  q0  and  n  are  assumed  to  be  ion  independent 
and  t  is  assumed  to  be  host  independent,  the  average  qQ  and  n  were  taken 
and  fixed.  The  process  was  then  repeated  with  minimization  with  respect 
to  t  only.  This  process  yielded  the  following: 

qQ  =  -1.09,  n  -  0.977  ,  (  191 ) 

and  the  T  values  were  well  fitted  by 

T  “  0.767  -  0.00896N  ,  (192) 

N 

where  N  is  the  number  of  f  electrons  in  the  configuration  4f  .  The 
values  of  on  used  in  the  above  were  not  varied  in  the  minimizing  process 
but  were  interpolated  from  the  values  given  in  equation  (185);  that  is, 

cr2  *  0.6902  -  0.0121N  , 

o4  =  0.09  (all  N)  ,  (193) 

O  «■  -0.04  (all  N)  . 
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TABLE  15.  PHENOMENOLOGICAL  Bnm  FOR  SIX  RARE-EARTH  IONS 
IN  CaW04  (all  in  cm-1) 


Ion 

820 

B40 

B44 

B60 

RB64 

IB64 

Nd 

549 

-942 

1005 

-17 

94  7 

1 

Tb 

468 

-825 

872 

-290 

595 

160 

Dy 

428 

-825 

978 

-7 

448 

2.5 

Ho 

436 

-664 

779 

-33 

558 

196 

Er 

404 

-685 

728 

12 

452 

164 

Trt 

417 

-689 

754 

17 

504 

359 

57 
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The  predicted  values  of  the  pnm(  ^  »q0fr|)  for  the  entire  rare- 
earth  series  are  given  in  table  16.  The  results  given  in  table  16  when 
compared  to  table  15  show  that  the  difference  between  the  derived 
Bnm( T 1 W ^  an<5  phenomenological  Bnm  is  greater  for  the  low-N  ions 

in  the  4fM  configuration.  This  may  be  a  defect  in  the  theory,  but  not 
enough  data  on  the  low-N  ions  are  available  for  analysis.  One  of  the 
significant  results  of  the  analysis  was  that  it  led  to  a  reanalysis  of 
the  spectrum  of  Tb:CaW04  with  a  different  interpretation  of  the  experi¬ 
mental  data  (Leavitt  et  al,  1974). 


TABLE  16.  DERIVED  CRYSTAL-FIELD  PARAMETERS,  B^T^on), 
FOR  4fN  CONFIGURATION  OF  TRIPLY  IONIZED  RARE-EARTH  IONS 
(all  in  cm-1  ) 


N 

Ion 

o 

C'J 

CO 

w 

o 

B44 

B60 

RB64 

IB64 

1 

Ce 

441 

-1429 

1462 

16 

1251 

52 

2 

Pr 

424 

-1224 

1253 

13 

996 

42 

3 

Nd 

408 

-1059 

1083 

11 

805 

34 

4 

Pm 

411 

-1017 

1041 

10 

764 

32 

5 

sm 

408 

-938 

960 

9 

676 

28 

6 

Eu 

408 

-887 

908 

8 

626 

26 

7 

Gd 

406 

-824 

843 

7 

559 

23 

8 

Tb 

424 

-856 

876 

8 

591 

25 

9 

Dy 

428 

-831 

851 

7 

563 

24 

10 

Ho 

418 

-756 

774 

6 

488 

20 

11 

Er 

415 

-707 

724 

6 

439 

18 

12 

Tm 

435 

-729 

746 

6 

454 

19 

13 

Yb 

434 

-701 

717 

6 

429 

18 

The  values  used  are 

02  -  0.6902  -  0.0121N,  Oj  »  0.09  (all  N), 
o6  »  -0.04  (all  N),  x  0.767  -  0.00896N , 
q  »  -1.09,  and  n  "  0.977. 


More  recent  work  on  CaW04  (Morrison  at  al,  1977)  obtains  the 
following  values  s 

0 2  =  0.6846  -  0.00895N, 

04  =  0.02356  +  0.00182N, 

T  =  0.75(1.0387  -  0.0129N), 

O  6  =•  -0.04238  +  0.00014N, 
q0  =*  -1.150,  and 
n  =  0.962. 


-  >■ 
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The  on  values  are  interpolated  from  the  calculations  of  Erdos  and  Kang 
(1972)  for  Pr^+  and  Tm^  +  .  The  factors  in  equation  (189)  were  combined 
so  that 

pn  =  <rn>(  1  -  on)/tn  ,  (195) 

and  the  along  with  the  x  are  given  in  table  17.  Thus  we  have 

aiJT»*O'T0  "  Pn*rm(el0'1')  *  (196) 


TABLE  17.  VALUES  FOR  T,  <rn>Hp  ,  on,  AMD  pp  FOR  4f" 
CONFIGURATION  OF  TRIPLY  IONIZED  RARE-EARTH  IONS* 


Ion 

N 

T 

<rZ>HF 

<lJ,>HF 

<r6>HF 

a2 

°4 

°6 

p2 

04 

06 

Co 

1 

0.7693 

0.3360 

0.2709 

0.4659 

0.6757 

0.0254 

-0.0422 

0.1041 

0.7536 

2.3417 

Pr 

2 

0.7597 

0.3041 

0.2213 

0.3459 

0.6667 

0.0272 

-0.0421 

0.  1756 

0.6464 

1.8754 

Nd 

3 

0.7500 

0.2803 

0.1882 

0.2715 

0.6578 

0.0290 

-0.0420 

0.1706 

0.5776 

1.5897 

Pm 

4 

0.7403 

0.2621 

0. 1655 

9.2247 

0.6488 

0.0308 

-0.0418 

0. 1679 

0.5339 

1.4213 

Sm 

5 

0.7306 

0.2472 

0.1488 

0.1929 

0.6398 

0.0327 

-0.0417 

0.1668 

0.5049 

1.3210 

Eu 

6 

0.7210 

0.2347 

0.1353 

0. 1686 

0.6309 

0.0345 

-0.0415 

0.  1666 

0.4836 

1.2503 

Gd 

7 

0.7113 

0.2232 

0.1237 

0.1477 

0.6220 

0.0363 

-0.0414 

0.1668 

0.4656 

1.1873 

Tb 

0 

0.7016 

0.2129 

0.1131 

0. 1287 

9.6130 

0.0381 

-0.0413 

0.1673 

0.4990 

1.1232 

Dy 

9 

0.6919 

0.2033 

0.1037 

0.1119 

0.6041 

0.0399 

-0.0411 

0.1681 

0.4341 

1.0614 

Ho 

10 

0.6023 

0. 1945 

0.0954 

0.0981 

0.5951 

0.0418 

-0.0410 

0. 1692 

0.4217 

1.0119 

Er 

11 

0.6726 

0.1865 

0.0883 

0.0074 

0.5861 

0.0436 

-0.0408 

0.1706 

0.4126 

0.9826 

Tm 

12 

0.6629 

0.1790 

0.0820 

0.0787 

0.5772 

0.0454 

-0.0407 

0.  1722 

0.4053 

0.9649 

Yb 

13 

0.6532 

0.1717 

0 .0753 

0.0681 

0.5683 

0.0472 

-0.0406 

0.1737 

0.3938 

0.9120 

*The  units  of  <rn>HF  and  pn  are  in  An. 


At  present  we  use  the  results  given  in  table  17  to  calculate  crystal- 
field  parameters  given  by  equation  (196)  and  use  these  parameters  as 
starting  values  to  best  fit  experimental  data.  We  generally  use 
Anmict'h)  with  n  =  1  in  the  process  (q  here  is  the  effective  charge  on 
the  ligands,  not  necessarily  oxygen).  After  obtaining  the  best-fit  Bnm, 
we  return  to  the  calculation  of  Anm(q)  and  vary  q  to  obtain  the  best  fit 
by  minimizing  the  quantity 

2  =  UBnm  “  Pn'W'*)]2  *  <197) 

nm 

Following  this,  we  obtain  the  Anm(q)  for  odd  n  and  use  them  to  calculate 
the  intensities  using  the  Judd-Ofelt  theory. 

At  present  we  have  not  included  the  dipole  contribution  to  the 
Anm(q)  but  intend  to  do  so  as  soon  as  possible.  The  old  n  in  the  three- 
parameter  theory  will  be  replaced  by  a,  the  polarizability  of  the  con¬ 
stituent  ions  in  low-symmetry  sites.  We  believe  that  this  latter  proce- 


dure  ( including  finding  new  p  values)  will  give  much  better  results 
than  obtained  by  the  older  theory.  In  our  projected  reanalysis  we  will 
have  the  good  phenomenological  Bnm  for  R^+sLaF.j,  R3+:LaCl3,  and 
R^+jLiYF^  (these  are  reported  by  Morrison  and  Leavitt,  1981),  and  will 
soon  have  R^  +  sY203,  in  addition  to  the  Bnm  for  R^+:CaWO^  used  in  the 
older  theory.  These  data  should  be  sufficient  to  form  a  very  stringent 
test  of  a  newer  three-parameter  theory. 


5.  CRYSTAL-FIELD  EFFECTS  NOT  YET  FULLY  INCORPORATED 


5.1  Self-Consistent  Point  Dipole  and  Point  Multipole 


In  section  4  we  discussed  the  point-charge  contribution  to  the 
multipolar  field  components  Anffl.  It  was  early  recognized  by  Hutchings 
and  Ray  (1963)  that  the  multipolar  components  of  the  constituent  ions 
contribute  to  the  A  at  the  si£e  occupied  by  the  rare-earth  ion.  For  a 
point  charge  eZ^  located  at  R^  from  a  rare-earth  site,  we  have  the 
electric  potential 

-*■  gZj 

♦(?)“7 *“T7  ‘  (198) 

I  | 


The  potential  energy  of  the  rare-earth  electron  at  r  is 
U(?)  »  -e<J>(r) 


(199) 


where  we  have  expanded  the  denominator  of  equation  (198)  in  the  spher¬ 
ical  tensors  discussed  in  section  2.  If  we  write  equation  (199)  as 


then 


U(?) 


1  A*„,rnC__(r )  , 

l  nm  run 


nm 


*(°) 

nm 


-«2 


l 


A 

z ,  c  (r,  1 

i  nrrr 


n+1 


(200) 


(201) 


where  the  sum  on  i  covers  all  the  ions  of  charge  eZ^  in  the  solid.  This 
result  we  derived  in  section  4,  expressed  in  slightly  different  form. 
It  seems  natural  to  extend  equation  (200)  to  the  form 
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(202) 


U<?)  ~  rl  *S>*r"Cn.(?>  ' 

nm 

k 

(k )  ± 

and  relate  the  A  m  to  the  various  k-pole  moments  of  ligands  at  R^.  To 

relate  the  the  multipole  moments,  eQ.  ,  we  need  first  to  express 

the  electric  potential  at  the  rare-earth  electron  due  to  the  multipole 

moment  eQ^q(i)  at  R^. 

The  electric  potential  due  to  a  multipole  distribution  at  is 


given  by 


$(?)  -  e  |  (-1  <r‘(m)k(q)  I n+k(m+q) > 


(203) 


x  Ea±JLtinta^?.i.)  rnc 

0n+k+ 1  r  cnm(r> 


where 


/  2k+2m  (  2k+2n)  I 

^  2n  >  "  (2n)l (2k) I 

(the  details  of  the  derivation  of  this  result  will  be  given  later). 
Thus,  since  U(£)  =»  -e<})(£)  we  find,  using  equation  (202),  that 

*  -e2  l  (_1  )1<CCa(i^  2°2nk^  ^  <n(m)k(q)  | n+k(m4q ) > 

q,i 


C  ,  f  R .  ] 
n+k,m+q'-  LJ 

n+k+1 


(204) 


If  we  let  k  *  0  in  equation  (204),  we  obtain  (if  k  -  0,  q  =  0) 


*42’  -  -2  l  &»>  ^4^  - 


(205) 


which  is  identical  to  the  result  given  in  equation  (201)  if  we  identify 
Qqq( i )  with  55^  (the  number  of  charges)  given  there.  The  result  for  k  = 
iis 

A*1*  -  e2  £  Q?  (2!J+2)1^2  <n(m)1(q)  |n+1  (m+q)> 
nm  L .  1qv  2n  '  ' 

q»i 


Cn+1  ,m+q^  Ri^ 


(206) 


since 


Qlq  =  (-1  > 

<n(m ) 1 (q ) j  n+1 (m+q  ) >  =  ( -1 ) 1 <1 ( _q )n+1 (m+g ) | n(m) > 

Then  with  these  substitutions  in  equation  (206),  we  have 

A(1)  =  -e2  V  ✓(n+1)(2n+37  Q,(  1 5  <  1  (q  )n+1  (m-q )  |  n(m) > 
q,  l 

A 

C  ,  C  R .  ) 
n+1  tin-q^  i; 


(207) 


which  is  indentical  to  the  result  previously  published  (Morrison,  1976), 
if  we  identify  eQ1  ( i )  =  P1q(i ) .  Thus,  if  we  knew  the  Q^.  ( i ) ,  we  could 
easily  calculate  the  by  using  equation  (204).  Unfortunately,  the 
real  difficulty  is  determining  the  Q^gd).  In  what  follows  we  shall 
restrict  our  discussion  to  the  dipole  case,  k  =  1,  and  let  =  pq  and 
express  the  results  in  Cartesian  vectors. 


At  sites  of  low  symmetry,  an  electric  field  can  exist  whose 
value  is  determined  by  the  various  point  charges  of  the  solid.  The 
electric  field  due  to  the  poinv  charges  of  the  solid  at  a  site  of  low 
symmetry  is  given  by 

_  q.R.  , 

,  -  l  -Ml  ,  ( 208) 

3  . 


R. 


and  the  field  generated  by  the  point  dipoles  is 


+d 


l 

i 


3R.  .IR.  , 
iJ 


(209) 


The  dipole  moment  at  site  j  is  then  given  by 


(210) 


where  a  is  the  polarizability  of  the  ion  at  site  j.  (If  more  than  one 
species  is  considered,  the  polarizability  of  each  type  must  be  used.) 
The  sum  in  equation  (208)  presents  no  problem  and  can  be  done  quite 
simply.  To  perform  the  sum  in  equation  (209),  it  is  convenient  to 
assume  a  fixed  coordinate  system  in  the  unit  cell  and  an  associated 
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reference  point  (say  position  1);  then  each  dipole  moment,  $y  can  be 
related  to  the  dipole  located  at  the  reference  moment,  p^ ,  byJthe  sym¬ 
metry  operation  of  the  crystal.  Similarly,  the  field  at  each  point, 
can  be  related  to  .  Having  done  this,  we  can  write 

,  (2ii) 

and  from  equation  (210) 


The  result  given  in  equation  (212)  can  then  be  solved  for  the  dipole 
moment  p^  to  give 

(213) 


The  result  obtained  in  equation  (213)  is  rather  interesting;  if  the 
polarizability,  a,  is  near  the  reciprocal  of  one  of  the  eigenvalues  of 
the  G  matrix,  then  the  dipole  moment  becomes  excessively  large.  This  is 
suggestive  of  the  type  of  catastrophe  that  occurs  in  the  onset  of  a 
ferroelectric  transition.  Such  a  situation  would,  perhaps,  be  modified 
by  the  inclusion  of  the  higher  multipole  moments  in  the  calculation.  It 
should  be  pointed  out  that  the  G  matrix  defined  in  equations  (209)  and 
(212)  is  dependent  only  on  the  lattice  constants  and  the  symmetry  of  the 
crystal.  The  results  here  were  expressed  in  terms  of  Cartesian  coordi¬ 
nates  but  can  equally  well  be  done  in  spherical  tensors.  If  higher 
moments  were  considered  the  spherical  tensor  form  is  much  more  conven¬ 
ient. 


*►  *>o 

pt  =  aB(1).E1 

where 

B  =  ( 1  -  oC)~ 
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TABLE  18.  SPACE  GROUP  88  (FIRST  SETTING): 
COORDINATES  OF  ALL  IONS  IN  A  UNIT  CELL  OF  YLiF4  AND 
DIPOLE  MOMENTS  OF  EACH  ION 
(px#  Py>  and  p„  of  site  1  are  chosen  as  u, 
v,  and  w  respectively) 


Site 

Ion 

X 

y 

Z 

Px 

py 

pz 

»kqa 

1 

F 

X 

y 

z 

u 

V 

w 

1 

2 

F 

y 

-X 

-Z 

V 

-u 

-w 

(-1)k(i)q 

3 

F 

-X 

-y 

Z 

-u 

-v 

w 

(-1)q 

4 

F 

-y 

X 

«2 

-V 

u 

-w 

<-1)k(-i)q 

5 

F 

1/2  +  x 

1/2  +  y 

1/2  + 

z 

u 

V 

w 

1 

6 

F 

1/2  +  y 

1/2  -  x 

1/2  - 

z 

V 

-u 

~w 

(-1)k(i)q 

7 

F 

1/2  -  x 

1/2  -  x 

1/2  - 

z 

-u 

“V 

w 

M)q 

8 

F 

1/2  -  y 

1/2  +  y 

1/2  + 

z 

-V 

u 

-w 

(-1)k(-i)q 

9 

F 

X 

1/2  +  y 

1/4  - 

z 

u 

V 

-w 

(-l)ktq 

10 

F 

y 

1/2  -  x 

1/4  + 

z 

V 

-u 

w 

<-i)q 

11 

F 

-X 

1/2  -  y 

1/4  - 

z 

-u 

-v 

~w 

<-1)k 

12 

F 

-y 

1/2  +  x 

1/4  + 

2 

-V 

u 

w 

(i)q 

13 

F 

1/2  +  x 

y 

3/4  - 

z 

u 

V 

-w 

(-1 )k+q 

14 

F 

1/2  +  y 

-X 

3/4  + 

z 

V 

-u 

w 

<-1)q 

15 

F 

1/2  -  x 

-y 

3/4  - 

z 

-u 

-v 

-w 

(-1)k 

16 

F 

1/2  -  y 

X 

3/4  + 

z 

-V 

u 

w 

(l)q 

17 

Li 

0 

0 

0 

- 

- 

- 

- 

18 

Li 

0 

1/2 

1/4 

- 

- 

- 

- 

19 

Li 

1/2 

1/2 

1/2 

- 

- 

- 

- 

20 

Li 

1/2 

0 

3/4 

- 

- 

- 

- 

21 

Y 

0 

0 

1/2 

- 

- 

- 

- 

22 

Y 

1/2 

0 

1/4 

- 

- 

- 

- 

23 

Y 

1/2 

1/2 

0 

- 

- 

- 

- 

24 

Y 

0 

0 

3/4 

— 

aTh«  last  column  relates  those  O^g  to  the  reference  point  0kg-  Thus 
the  Ojt^  for  fluorine  are  all  relatedrto  site  1. 


5.2  Self-Consistent  Results  for  Scheelite  Structure 


The  procedure  given  above  is  rather  involved,  so  we  shall  go 
into  the  derivation  of  the  G  tensor  for  the  Scheelite  structure  (CaW04, 
LiYF4,  etc).  The  space  group  for  Scheelite  is  88  in  the  International 
Tables?  the  position  of  all  the  constituents  is  given  in  table  18.  To 
be  specific,  LiYF4  has  been  chosen,  and  the  fluorine  in  site  1  at  x,  y, 
and  z  has  been  chosen  as  the  reference  point  for  the  dipoles  (u,  v,  w) 
and  all  other  dipoles  in  the  unit  cell  are  related  to  u,  v,  and  w.  No 
dipoles  can  oxist  at  the  Y  or  Li  sites  since  the  lowest  fields  at  these 
sites  are  quadrupole  (k  -  2). 


MBKiirMiTilTITKiBBW 


To  evaluate  G  for  the  Scheelite  structure,  we  shall  choose  the 
ion  at  site  1  in  table  18  as  j  in  equation  (209).  The  ft  in  equation 
(202),  including  the  translational  vectors  (£  in  x,  m  in  y,  n  in  z),  is 


.  /.A 

®i ,  i  =  { S.+x±--x)aex  +  (m+y^-y)aey 


(214) 


+  (n+Zj-z)ce„ 


where  we  shall,  during  this  discussion,  suppress  the  explicit  dependence 
of  R  on  l,  m,  and  n.  We  write  equation  (209)  as 


+d  -t-d  ,  >d 

E,  “  Fi  + 


(215) 


where 


+d  r- 

F  =  -  ) 


-► 

P4 


i  R 


i,1 


(216) 


and 


Fd - a- 

x  „3 


u 


Then  we  can  write 


(-1) 


X  4  <  T.3 

i“1  R_ , 


*d 

ri  =  3 

ft.  •  p. 

r  ♦  ii  i 

(217) 

f  Ri,i 

i 

t 

,  m,  and  r. 
f  as 

are  implicit. 

Then  using 

table 

18,  we 

v  +  — 

+  v  _ 

u 

v  u 

_3  3 

R2 , 1  R3 , 1 

R4,1 

RS.1 

R3  R3 

6,1  7,1 

u  V 

i  u 

+  V 

u 

V 

R3  R3 

R9,1  R10, 

i  Ri  i ,  i 

R3 

12, 

R3 

1  R13,1 

R3 

14,1 

V 

R3 

R16,1 

(218) 

i  8 

_ „  .  r 

( 0  )w 

(219) 

—  u  +  i 
M  i'1 

R3 

2i,1 

• 
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If  we  let  G 


and  restore  the  %,  m,  n  sum,  we  have 


-  G'  +  3G" 


A* 


G'  = 
XX 


G'  = 

xy 


8 

l  l 


8 

l  l 


(-1)' 


£,m,n  i=1  R 

(-1) 


2i+1 , 1 
i 


£,m,n  i=1  R' 


21 , 1 


G*  =0  . 

xz 


(220) 


(221) 

(222) 


By  similar  methods  we  obtain 

r  r  i  1  1 

G’  =  l  1  (-D  --3 -  »  (223) 

i,m,n  i=1  Ri  Ri+8 

and  the  G’  tensor  is  symmetrical. 

To  evaluate  the  G"  term,  the  procedure  is  precisely  the  same  as 
to  evaluate  the  G1  term,  except  that  we  relate  this  to  equation  (217). 
It  is  convenient  to  express  R^'P.  explicitly  in  tabular  form  as  given 
in  table  19  for  easy  reference  when  writing  out  each  term  of  G" .  We 
shall  not  write  out  the  detailed  expression  as  in  equation  (218),  but 
this  procedure  gives 


G" 

xx 


G 

xy 


G" 

xz 


G" 

yz 


^i+1  ”  X2±Y2i^  ' 

J  <-,)1[x|i  -  . 

'  W».l  ' 

-J,  -  Wi+sl  • 


(224) 


(225) 


(226) 


(227) 
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(228) 


z*  -  -j,  -  z!+3I  • 


where 


x  y  _  a  x(i)y(i) 

1  1  R5 

i,1 

2  2 
„2  a  x(i >  Z 

X.  --5  . 

i,1 

„  2  ,  .  .2 

Y2  -  a  iLll 

L  ^  1  ' 
1/1 


v  _  acx ( i ) z ( i ) 

Vi  ■  5 

i/1 


and  all  of  the  sums  in  equations  (224)  through  (228)  have  the  sum  over 
l,  m,  and  n  implied.  The  G"  is  symmetric  (this  can  be  shown  directly 
from  evaluating,  for  example,  G£y  and  G^  independently). 

TABLE  19.  VALUES  OF  R*P  FOR  DIFFERENT  SIRES  IN  SCHEELITE3 


Site 

Px 

py 

Pz 

R  •  P 

1 

u 

V 

w 

x<  1  )u  +  y( 1 )v  + 

z(  1  )w 

2 

V 

-u 

-w 

x( 2 )v  -  y(2)u  - 

z(2)w 

3 

-u 

-V 

w 

-x(3)u  -  y(3)v  + 

z(  3)w 

4 

-V 

u 

-w 

-x(4)v  +  y(4)u  - 

z(  4  )w 

5 

u 

V 

w 

x(5)u  +  y(5)v  ♦ 

z  ( 5  )w 

6 

V 

-u 

-w 

x(6)v  -  y(6)u  - 

z(6)w 

7 

u 

-V 

w 

-x(7)u  -  y(7)v  + 

z(7)w 

8 

-V 

u 

-w 

-x(8)v  +  y(8)u  - 

z(8)w 

9 

u 

V 

-w 

x(9)u  ♦  y(9)v  - 

z(9)w 

10 

V 

-u 

w 

x( 10)v  -  y( 10)u 

+  z ( 10  )w 

11 

-u 

-V 

-w 

-x(11)u  -  y( 1 1 )v 

-  z( 1 1 )w 

12 

-V 

u 

w 

-x(12)v  +  y ( 12 )u 

-  z(12)w 

13 

u 

V 

-w 

x( 1 3 )u  +  y(13)v 

-  z( 13)w 

14 

V 

-u 

VI 

x{ 14)v  -  y ( 14)u 

+  z(14)w 

15 

-u 

-v 

-w 

-x( 15)u  -  y( 15)v 

-  z( 15)w 

16 

-V 

u 

w 

-x(16)v  -  y(16)u 

-  z ( 1 6)w 

*x(l) 

-  t  +  s 

i  "  *V 

- 

m  *  yi  -  yir  z (i) 

»  n  +  -  Zj 
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The  equations  for  G'  and  G"  were  Calculated  for  several  lat¬ 
tices  and  the  results  are  given  in  table  20.  The  crystal^ axial  field 
components,  Ano,  were  computed  for  ^aW04  using  a  -  2.4  A'  and  oxygen 
charge  of  -2e,  and  using  a  =  0.24  A  and  oxygen  charge  of  -e.  The 
results  are  shown  in  table  21  (Morrison,  1976). 

d  After  the  above  work  had  been  done,  the  dipole  terms  in  the 

A  were  programmed  for  a  computer  for  all  the  230  space  groups.  In 
the  program  any  number  of  inequivalent  sites  can  have  an  associated 
dipole  moment  (we  only  considered  one  type  of  site  above).  Recently  the 
members  of  P.  Caro's  group  in  France  and  G,  F.  De  S&'s  group  in  Brazil 
(Faucher  and  Malta,  1981)  have  included  the  dipole  and  quadrupole  moment 
in  a  self-consistent  manner  for  LaCl3  and  have  found  that  with  the 
reported  values  of  the  dipole  and  quadgupole  polarizabilities  the  re¬ 
sultant  Aq  is  much  larger  than  A°  or  A 

nm  nm  nm 

5.3  Self-Induced  Effects 


When  a  rare-earth  ion  i3  immersed  in  a  solid  it  is  possible  for 
its  electrons  to  experience  a  field  due  to  the  reaction  of  the  medium 
back  on  the  electrons.  Both  this  type  of  field  and  the  external  fields 
due  to  the  point  charges  of  the  medium  can  exist.  This  reaction  is 
identical  to  the  classical  problem  of  a  charged  particle  interacting 
with  its  induced  image  in  a  conducting  plate  or  sphere.  The  interaction 
was  recognized  by  Judd  (  1977),  and  it  was  he  who  suggested  the  polar¬ 
ization  of  the  ligands  as  «.  possible  source  of  a  two-electron  crystal- 
field  interaction.  In  this  section  we  shall  consider  the  development  of 
this  interaction  as  derived  earlier  (Morrison,  1980),  using  the  same 
technique  used  in  the  earlier  work.  In  later  sections  this  interaction 
will  be  developed  in  a  more  general  way,  deriving  the  multipolar  inter¬ 
action. 

We  consider  an  electron  at  r  on  a  rare-earth  ion  and  a  ligand 
at  R  with  polarizability  a.  The  electric  potential  created  by  the 
electron  is 


♦  = 


I  R-r  | 


(229) 


The  electric  field  at  the  ligand  is 


E  “  "VR*  , 


where  VR  indicates 
to  R.  Than, 


that  the  derivative  should  be  taken  with  respect 


E  = 


,  +  •* 
-a (R-r) 

i  A  1 3 
R-r 


(230) 
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TABLE  20.  G  TENSOR  AND  X-RAY  DATA  FOR  SEVERAL  COMPOUNDS  (1  A3  ) 


w 

u 

5 

d 

> 

o  u 

s  s 

* £ 
o  o 


j  a 
M  X 
W  X 

Du  O 

0  ° 

sS  SSr 


Sr 

rn  »  . 


o 

X 

X 

o 

bu 

o 

x 

fr* 


8  S 

d  s 

3  ^ 
3  2 

+  tu 
t-  o 
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►J 


s 


X 

o 

a4 


*9 

<N 

o 

ft 

0 


0) 

(N 

I 


X 

o 

O4 


*< 

fN 

I 

0 


O 
•d  n 


o 

v 


o 

T3  C 
4. 


°*c 


4J 

G 

0} 

c 


o 

o 
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PMAW 


The  dipole  moment  induced  on  the  ligand  is  given  by 

P  =  a&  ,  (231) 

where  a  is  the  polarizability  of  the  ligand. 


Now  if  we 

electric  potential  at  point  R  from  that  origin  is 


consider  a  dipole  at  some  arbitrary  origin,  the 


p*R, 


R 


To  find  this  potential  at  the  electron  itself,  we  let  R 
equation  ( 232 )  becomes 

*  =  7P'<x-n  . 

1  +  i3 

R-r 


=  -(R-r) 


(232) 


Then 


(233) 


The  energy  of  the  electron  interacting  with  this  potential  is  given  by 

.+  -e 


U(r,R)  =  —  <j>1  ( r ) 

-*•  ,+  +  . 
e  p « ( R-r ) 

2  i&  -*-|3 
R-r 


(234) 


where  the  1/2  is  due  to  a  self-interaction.  We  can  write 


R"r  =  -v  1 


I R-r I  3 


R  +  i 

I  R-r  I 


(235) 


Then  equation  (234)  becomes 


■*  ■+■  -e  > 

U(r,R)  =  —  p  •  VR 


I  r  T  , 

I  R-r  | 


(236) 


and  similarly 


E  =  eV, 


B  '5-fi 


(237) 


Using  the  result  of  equation  (237)  in  equation  (231)  and  substituting 
the  result  into  equation  (236),  we  have 

2  / 


1 


U(?,5)  =  -f-  «JVR 

\  I R-r | , 


1 


R  i  +  1 

I  R-r  1/ 


(238) 
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where  VR  operates  only  on  the  function  on  its  immediate  right.  To 
further  reduce  the  result  given  in  equation  (238),  we  consider  the 
operation 

?2(V2)  -  ^V2>2  +  <|»272^  +  2(V*1MV*2)  •  (239) 

If  ^  and  i)»2  satisfy  Laplace's  equation  (which  they  do),  then 

^(*,*2)  a  2(7*1)-(Vi|;2)  .  (240) 


If  we  identify^ij^  and  ^  with  1/|R  -  r|  in  equations  (238)  and  (240),  we 
can  write  U(r,R)  as 

0(L*)  -  2Hi  vj .  (241) 

To  proceed  further  we  must  expand  the  factors  on  the  right  side 
of  equation  (241).  First  we  notice  that 


If  we  let 


then 


A  A 

with  z  =  r»R. 


,♦  +1 2 

R-r 


2  2  _+ 

R  +  r  -  2r*R 


-  2RrKsr  - • 


t  = 


2  2 
R  +r 

2rR 


t  +,2  2rR  "t-z ^ 


R-r 


The  expansion 


th  =  l  (  2  n+ 1 )  Qn  ( t )  Pn  ( z ) 


is  given  by  Rainviile  (I960);  the  leading  term  for  large  t  is 


„  .  2n(n! )2 

Qn(t)  “  -7TT - L - 

n  .n+i , _  .... 

t  ( 2n+1 ) I 


(242) 


(243) 


(244) 


(245) 


(246) 


From  equation  ( 243 )  we  have 

_2n+1,  ..2  n+1 

«...  2  (n! )  r 

0  (t)  “  . . . .  —  — > 

wn'  ;  ( 2n+1 ) I  n+1 

R 

for  large  R. 

Substituting  the  result  of  equation  (247)  into  equation 


gives 


_2n.  , . 2  n 

i  v>  • 


s  R-r  r  ( 2n ) ! R 

From  the  Legendre  addition  theorem,  we  have 

v*)  - 


and 


nm 

,2n,_,v2  _n 


I  Cnm<*>CZmW 


| R-r | 2  nfm  (2n)1  Rn+2  Cnm<r>Cnm<R> 

2 

The  remaining  necessary  operation  is  V  ,  which  can  be  written 

R 

v2  =  1_  d_  f  2_d^  _  (L)l 

k  2  dR  dR'  2 
R  R 

The  only  term  in  equation  ( 250 )  that  this  operates  on  is 

^n+2  Cnm(R)  * 

Using  equations  (251)  and  (252),  we  have 

V2  fr-n-2c*  (R)  1  =  R l-l.Lnj‘2 )  _  n  ( n+1 )  ~[  p*  ,  pv 
VR  Lr  Suit  k'  J  [_  Rn+4  Rn+4  J  LnmlR) 


Finally, 


72rR”n,“2C*  (R)  1  =  2. ( U+.l  )  r;*  (d) 
vRlR  cnmtR'-i  Rn+4  LnatR'  ' 


where  we  have  used 


(£)2C  (R)  -  n(n+1)C  (R)  , 

nm  nm 


(247) 

(244) 

( 248 ) 

(249) 

(250) 

(251) 

(252) 

(253) 

(254) 

(255) 
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a  result  we  discussed  In  section  2.1.  The  result  in  equation  (254) 
substituted  into  equation  (250)  gives 


72— J--  I 

R  .  ,2  L 


l  _4  c  ,;)c. . 

L  (2n)!  n+4  nm  nm 

nm  R 


(256) 


The  result  given  in  equation  (256)  is  substituted  into  equation  (241). 
This  result,  when  summed  over  all  ligands  at  with  polarizability  ay 
produces 


_2n+1  , .  , , , 
2  nl(n+1) 

( 2n) ! 


1  n 

J  nmv  3'  n_  ,  4 
~  ~T+4  r  Cnm  r  * 
R1  J 


(257) 


If  we  write  equation  (257)  as  we  have  previously  done  with  the  point- 
charge  model. 


u(r)  =  J  (w  rIlcWr>  ' 


(258) 


we  have 


»S1  -  (4) 

nm  \  4/ 


22n+1nl  (n-M)  I 
(2n)  ! 


a  ,C  (r.1 

nm'-  -W 


(259) 


which  are  the  self-induced  crystal-field  components  due  to  induced 
dipoles  only.  Higher  order  multipole  moments  can  be  induced  on  the 
ligands,  and  these  multipoles  will  contribute  a  correction.  From  pre¬ 
vious  experience,  we  should  anticipate  the  total  self-induced  multipole 
fields  to  be  of  the  form 


l  AS1(k) 
k=*1  nm 


(260) 


with  the  result  above  being  A  ( 1 ) . 

nm 


As  in  the  point-charge  model,  if  we  express  all  lengths  in  A 
and  Oj  in  A3  then  equation  (259)  becomes 


.SI,,,  o  8 

A  ( 1 )  -  -  —  x  10 
nm  \8ir 


„2n+1  , ,  . ,  a  .C  (r , 1 
2  nl (n+1 ) 1  r  j  nm' 

(2n)i  5  _n+4 

j  Rj 


(261) 


to  express  A  in  units  of  cm“VAn  (o0/8tt  x  10®  =  29,035). 


73 


6.  MISCELLANEOUS  CRYSTAL-FIELD  EFFECTS 


6.1  Judd's  Interaction  for  Two  Electrons 


The  interaction  considered  hare  is  a  development  of  a  sugges¬ 
tion  by  Judd  (1977)  concerning  a  possible  origin  of  two-electron  crys¬ 
tal-field  effects.  Specifically,  Judd  suggested  that  such  terms  would 
arise  if  one  of  the  electrons  in  the  configuration  4fN  polarized  a 
nearby  ion,  and  the  remaining  N-1  electrons  interacted  with  the  induced 
multipolar  moments.  The  investigation  of  this  interaction  was  performed 
later  (Morrison,  1980),  assuming  only  a  dipole  polarizability.  The 
interaction  for  two  electrons  that  resulted  is 


V(  1 , 2 , R)  -  l  F(abk)r*Caa(;i)r5cbfq_a(r2) 

€1  Ct 

b,k,q 

A 

C*  ( R ) 

X  <a( a)b(q-a) |k(q) > 

R 

where 

F(abk)  =  “  §-^<a(0)b(0) |k(0)> ( (a+b+1 ) (a+b+2)  -  k(k+1)] 

and  a  is  the  dipole  polarizability  of  the  ion  at  R. 


(262) 


9 


The  development  of  the  result  given  in  equation  (262)  was 
similiar  to  that  given  in  the  derivation  of  the  self-induced  field  in 
section  5.3.  For  the  full  multipolar  result  we  shall  use  more  general 
methods . 


The  electric  potential  of  an  electron  at  r,  as  seen  at  a  ligand 


at  R,  is 


where 


♦tig -7? 


-e 

R.  I 


R,  -  R  -  r, 


The  multipolar  inducing  field  Enm  at  R^  can  be  defined  by 


<j>fR.+x)  =  -  T  E*  xnc  (x) 
1  1  L  nm  nm 

nm 

By  expanding  equation  (263),  we  obtain 


e  ^  (.^n  £nmWl_  Cnm(x)  , 

nm 


'1 


(263) 


(264) 


(265) 
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then  comparing  equation  (265)  with  equation  (264)  gives 

E  _  -  e ( - 1 ) nC  (r.  )/r!J+1  • 

nm  nm'-  1 '  1 

The  multipole  moment,  Q  _,  is  given  by 


Q__  =  a  E  , 
nm  n  nm 


(266) 


(267) 


where  the  multipole  polarizability  is  a 


n 


The  electric  potential  at  an  arbitrary  point  R3  from  a  multi¬ 
pole  distribution  is  given  by 


Q *  C  (r  ) 
+  ->  p  nm  nm\  3/ 


♦B(V  -  I 


n+1 


(268) 


and  the  energy  of  an  electron  at  r^  interacting  with  the  multipoles  is 


0  .  -e*Q  (J3  -  -S2) 

with  r  -  $  -  ?2.  From  equation  (268),  we  obtain 


r  QSm(“1)nC„m(R0l 
nm  R- 


(269) 


(270) 


Now  from  equations  (267)  and  (266)  we  have 


0  =  ea_ (-1  )n 

wnm  ectn'  1;  Dn+1  ' 


(271) 


which,  when  substituted  into  equation  (270),  gives 
»(*,M)  ■  -2  L*n  • 


nm  “  R  J 


(272) 


If  we  write 


we  have 


,r<(n)  (■*■  *  +  -\  t  p 

U  (rtr2,R)  =  ~®  an  I 


.  c*  (r, )  c  (r.) 

2  p  nm'1  1'  nm'-  2' 


m  R 


n+1 


n+1 


(273) 


(274) 
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•A  >. 


If  we  were  considering  the  self-interaction,  at  this  point  we  would 
let  $2  =  and  take  half  the  results.  The  sum  on  m  would  then  collapse 
to  unity  and  U( n) (r1 ,r1 , r)  -  -e2(  0,^/2  )/R2n+2. 

However,  the  two-electron  interaction  is  more  complicated.  We 
use  the  two-center  expansions  to  obtain 


Cnm^1  j  _  v  f2a+2 n'^2 


1Caa(r1 ) 


-  L  (  2an)  <a(a)n(m)  |a+n(a+m)>  C»+n>a+m(R) 


aa 


(275) 


where  R1  =»  R-r  ^ ,  and 
Cnm(^2)  _  r  (2b+2n)l/2 

Rn+!  -  ^  l  2b  J 


<n ( m ) b (  3 )  I b+n (  3+m )  >  cb+n,3+m(R) 


(276) 


where 


-►  >  + 
R2  "  R-r2 


As  indicated  in  equation  (272),  equations  (275)  and  (276)  are 
to  be  multiplied  together.  When  these  two  equations  are  multiplied,  the 
two  spherical  tensors  in  R  can  be  recoupled  as 


C*  (R)C  (R)  =  ( - 1 ) °t+m  l  <b+n (0)a+n(0)|k(0)> 

a+n ,  ot+m  b+n ,  3+m  £ 

A 

x  <b+n( 3+m)a+n(-a-m) |k( 3-o) >C,  (R)  , 

k,  3-a 


(277) 


where  we  have  used 


C* 

a+n,  cc+m 


A 

(R) 


( - i ) a+mc 


a+n, -a-m 


(R) 


(278) 


It  should  be  noted  that  the  resultant  projection  in  equation  (277), 
lCk  6-a(R*  J'  *-s  ^dependent  m*  Thus  with  a  proper  recoupling  of  the 
C-G  in  equations  (275)  and  (276),  the  sum  over  m  can  be  performed. 
Selecting  the  independent  terms  from  the  product  of  equations  (275)  and 
(276)  and  the  result  of  equation  (277),  we  have 


S  -  l  (-1 )m+a<a( a)n(m) |a+n( a+m) ><n(m)b( 3) |b+n( 3+m)> 

m  (279) 

x  <b+n( 3+m) a+n (-a-m) | k( 3-a) >  » 
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which,  when  further  reduced,  gives 

U<">(f,f2,a)  .  -e2o„  Io|  <b*„(0).+n(0)|k(0)>[(2»*2")(“J2“)]1/2 

b0  (280) 


x  r 


aC  (r  )rbCt  ^R2^SCk,B-a(R) 
1caa>-r1  Jr2cb0 „a+b+2n+2 


Thus  the  final  desired  result  is  obtained  if  we  know  S. 
(279)  we  rearrange  the  C-G  as  follows: 


In  equation 


<a(  ct)n(m)  |a+n(a+m)> 
<n(m)b( 0) | b+n( 0+m)> 


( -1 )a"a  (2Z±2n+J ]l/2  <a(-a)a+n(  ort-m) |n(m)>  , 

1/5  (281) 

(-!)”  <n(m)b+n(-0-m)|b(-0)>  . 


We  then  recouple  the  two  C-G  on  the  right  to  give 


<a(-a)a+n(or*-in)  |n(m'Xn(m)b+n{-0-m)  |b(-0)> 

"  I  / ( 2f+1 ) ( 2n+1 )  W(a,a+n,b,b+n;nf )<a+n( a+m ) b+n ( - 0-m ) | f ( a-0) > 


x  <a(-a)f (a-0) |b(-0)>  , 

The  sum  on  m  can  now  be  performed  (note  that  the  phase,  (-1)m,  in  eq 
(281)  cancels  the  ( — 1 )m  in  eq  (279))  if  we  change  the  phase  in  the  first 
C-G  on  the  right  side  of  equation  (282).  This  then  fixes  the  sums  on  f 
at  k.  Thus 


S  *  (-i) k~b+n [ I3&±3gi{ 1 2b*2n+ 1 1] 1  /2  /-^FTr 

x  W(a,a+n,b,b+n;nk)<a(-a)k(a-0) |b(-0)> 


(283) 


The  C-G  in  equation  (283)  can  be  rearranged  to  give 

l/2 

S  -  (-1)0l[(2a+2n+1)(2b+2n+1)  ]  W(a,a+n,b,b+n»nk)<a(a)b(  0)  |k(a-0)>  . 

(284) 
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If  we  let 


Pn(abk)  =  - ( c^e2 ) <a+n( 0 )b+n( 0 ) | k ( 0 ) >/ ( 2a+2n+ 1 ) ( 2b+2n+1 ) 

x  W(a,a+n,b,b+n;nk)[(2aJ2n)(2bj2n)]l/2  , 


(285) 


then,  substituting  into  equation  (280),  we  have 
,(n)  f±  ±  a  b 


U  (r1r2,R)  =  l  Fn<abk)r1r2  J  <a(oi)b(q-a)|k(q)>Caa(r1)cb  (r2) 

a,b  ~  ' 

k,q 


a 


x  cka|RI 
X  q+b+2n+2 


(286) 


which  is  the  final  form  of  the  two-electron  multipolar  interaction.  To 
obtain  the  result  given  in  equation  (262),  we  would  have  to  relate 
<a+1(0)b+1(0) |k(0)>  to  <a(0)b(0) |k(0)>  and  evaluate  W(a,a+1 ,b,b+1 , Ik) , 
both  of  which  can  be  found  in  Rose  (1957s  47,  227).  If  this  is  done, 
then  equation  (286)  will  reduce  to  equation  (262).  In  a  solid  the 
ligands  at  R  are  such  that,  when  the  sum  is  performed  over  the  ligands, 
only  certain  k  and  q  survive.  Much  of  the  above  derivation  has  been 
given  by  Judd  (1976)  in  a  different  context,  and  many  of  his  elegant 
techniques  could  be  used  to  simplify  the  resulting  expressions.  For 
example,  using  Judd's  notation  (1975),  equation  (286)  becomes 


-  l’ VabklrVXK)Sb(r2)WR,/R’ 

ab 
k 

where 


a+b+2n+2 


(287) 


[£a(rl)£b(r2)]kq  “  l  <a<a)b(q-a)|k(q>>Cao(r1)cbfq..0(r2)  • 

The  tensor  in  orbital  space,  given  in  equation  (286), 


Tka(a,b)  ”  ^  <a(a)b(q-a)  |k(q)>C  (r^t^  (r2)  '  (288) 
4  a  '4 

should  be  considered  carefully.  For  a  fixed  value  of  k  the  number  of 
terms  in  the  sum  over  a  and  b  is  restricted  by  a  +  b  <  k;  for  equivalent 
electrons  a  and  b  are  restricted  to  even  integers;  and  for  0  <  (a,b) 
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<  6,  the  total  number  of  terms  is  not  excessive.  But  since  a  and  b  can 
reach  the  maximum  value  of  6  for  the  configuration  4fn,  the  value  of  k 
in  the  k  sum  (similiar  to  the  lattice  sum)  must  go  up  to  12,  that  is,  k 
i  12. 


If  as  in  previous  work  (Morrison,  1980)  the  sum  over  all  the 
electrons  is  performed  in  equation  (286)  along  with  the  sum  over  the 
ligands,  the  results  are 

H(n)(R)  -  1/2  l  U(n)(r. ,r  ,R)  ,  (289) 

ij  3 

R 

where  the  factor  1/2  accounts  for  the  self-interaction  terms  that  are 
present  when  an  electron  interacts  with  its  own  induced  multipole,  as 
well  as  for  the  interactions  that  occur  twice  when  i  *  j.  This  inter¬ 
action  contains  a  large  number  of  corrections  to  the  free-ion  param¬ 
eters,  a  few  of  which  shall  be  discussed  in  the  following. 

6.2  Slater  Integral  Shifts 

The  Slater  integrals  for  the  free-ion  interactions  are  given  by 
the  Coulomb  interaction  as 


H  .  I  —2 - , 

i>j  ki-rj| 


which  for  equivalent  electrons  can  be  written 

,(k) 


where 


,<k> 


k,q 


°2  rt  rB  -u- 

r  > 


(290) 


(291) 


Since  the  interaction  represented  by  equation  (291)  is  spheri¬ 
cally  symmetric  in  the  space  of  all  the  electrons,  corrections  to  the 
can  only  arise  from  terms  in  an  interaction  that  are  spherically 
symmetric  in  the  space  of  the  electrons.  Thus,  in  equation  (286)  if  we 
let  k  *  0,  we  have  such  an  interaction,  and  the  following  results  are 
achieved: 
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(292) 


FnUbO)  =  -ane2<a+n(0)a+n(0)  |0(0)>/(2a+2n+1 )  (2b+2n+1) 
x  W(a,a4n,a,a+n;nO)  [ (2^2n ) (2b^n)  ] ^ 

=*  -one2(-1)a+n(2|^2n)//2a+2n+1  , 


where  we  have  used 

<a(a)b(-a) |0(0)>  =  (-1 )a“'V/2b+1  6  . 

an 


(293) 


and 


W(a,a+n,a,a+n;nO) 


( — 1 ) n/ 1 ( 2a+ 1 ) ( 2a+2n+ 1 ) ] 


1/2 


Then  equation  ( 286 )  becomes 

»<">(?,, f2.*)  -  -v2 i  (2*:2n)  i  'M 

a  a 

which  is  the  same  form  as  equation  (291). 


p,2a+2n+2 

Thus, 


(294) 


AF 


(k) 


-  'I 

i 


a  (i)e 
n 


('2a+2n-»  [<ra>] 2 


2a 


R2a+2n+2 

i 


(295) 


for  the  Slater  integral  shifts  due  to  the  electron  multipolar  interac¬ 
tion  with  the  ligands  of  multipolar  polarizabilities  a  . 

n 

6. 3  Shifts  of  Energy  Gap  4fN-4fN~1nft 

The  shifts  in  the  energy  gaps  between  the  4fN-4fN-1n£  dis¬ 
cussed  here  are  dependent  on  one-electron  operators,  and  it  is  best  to 
begin  the  discussion  by  returning  to  equation  (274)  and  proceeding  from 
there.  As  mentioned  following  equation  (274),  the  self-interaction  is 
obtained  by  letting  =  li  in  equation  (274),  obtaining 


U(n)U,R) 


e2a 


„2n+2 


(296) 


where  ■  £  -  r.  The  quantity  R^nf  in  equation  (296)  can  be  expanded 
as  in  equation  (242): 


R2n+2  =  £R2  +  j-2  _  2ft*?]n+1 


(297) 
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or 


R2n+2  =  (2Rr)n+1[t-z]n+1 


(298) 


where 


and 


t  =  (r2  +  r2)/2rR 


z  =  r*R 


Then  we  have 

.  ,  e  a  . 

ulnl(f,S)  -  -  "  1 


the  expansion 


2  __  n+1  .  .n+1 

( 2Rr )  |  t-z  | 


ti*  l  I2K*’>fik(t)Pk(2) 


is  given  in  Rainville  (1960:  182).  Also  we  can  write 
1 


(-1 )n  dn  1 


,n+1  n!  ,.n  t-z 
(t-z)  dt 

For  our  purposes  here  we  want  only  the  k  =  0  term  of  equation 
(P0(z)  =  1).  From  equation  (301)  we  know 

2 

„<n>  {r  _  ~e  an  .  1  <-Dn  Q(n)  . 

<r'R)  2  n+1  nl  C0  (t)  ’ 

( 2rR) 

By  using  equation  (303)  in  equation  (302)  with  k  =  0,  we  get 

Q0(t)  =  t”1  +  t'3  +  .  .  . 

and 

Qjn)(t)  *  (-1 )nnlt”n_1  +  (-1)n  (5t|)L  t“n"3  +  .  .  . 

Using  equation  (306)  in  equation  (304),  we  have 


2 

.  .  -e  a 

U  n  (r,R)  - 


2  ( 2rR)n+1  Lt 


1 _  (n-M)(n+2) 

n+1  n+3 

Ot 
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(299) 

(  300) 
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(302) 

(  303) 

(  302) 

(304) 

(305) 

(306) 

(307) 
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For  large  R  >  r,  equation  (299)  gives 


Thus  equation  (307)  becomes 


2  r  o' 

u<n>(r  R)  _  IfJil  __L_  +  _L_ 

U  (r,R)  2  R2n+2  3  R2n+4 


(  308) 


The  result  given  in  equation  (308)  is  a  generalization  of  the 
result  (Morrison,  1980)  for  a  ligand  with  dipole  polarizability  (n  =  1 
in  equation  (308)).  Thus  we  have,  for  n  =  1, 


.ala  4  ♦2gi 

R '  R 


(  309) 


where  we  have  written  =  a.  In  the  remainder  of  this  section  we  shall 
discuss  only  the  dipole  part  of  the  interaction  as  given  in  equation 
(309). 

If  we  let  w4jf  represent  the  lowest  energy  of  the  4fN  configu¬ 
ration  and  let  wni^i  represent  the  corresponding  energy  of  the  4fn_  n ‘ 
configuration,  we  may  write  W4f  =  V1°4f  +  A4f  and  ^'Z'  *  <  +  An'Z" 

where  and  are  the  corresponding  quantities  in  the  free  ion 

and  A and  A  ...  are  the  shifts  represented  by  equation  (309).  Thus, 
4f  n  Z  ' 

we  can  write 


(2  2  \ 

(^r  >f  0  I  ~  <r  *  A f  )  ' 


(  3  10) 


where  <r2>f  =  < 4fN_1n - Z ’ | r 2 | 4fN_1 n 'Z 1 >  and  <r2>4f  =  <4fN | r2 | 4fN>. 
The  result  in  equation  (310)  is  for  a  single  ligand  at  R,  and  if  we  let 
there  be  Zi  ligands  it  with  polarizability  o^,  we  have 


aiV 


(311  ) 
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where  a2  ••  <r2>U"  -  <r2>4f.  Since,  in  general,  <r2>f^i  is  greater  than 
<r2>4£,  we  have  the  result,  observed  experimentally,  that  the  free  ion 
interconf igurational  energy  gap  is  reduced  when  the  ion  is  embedded  in  a 
solid.  This  result  is  intuitively  obvious,  since  the  presence  of  the 
polarizable  ligand  tends  to  counteract  the  central  potential  of  the 
rare-earth  ion. 


Recently,  Yang  and  DeLuca  (1976)  have  measured  the  energy 


difference  of  the  lowest  4fN  and  4fN"15d  configurations  for  Nd3+, 


Er3  , 


LuP,,  and  LiYF.  (table  22), 


and  Tm3+  in  the  host,  materials  LaF^,  YF^,  uu,  3 
The  distances  of  the  nearest  neighbor  ligands  are  given  in’ table  23  for 
the  compounds  reported  by  Yang  and  DeLuca.  We  computed  the  values 
of  S  *  £jpiZ^e2/Rf  from  the  results  of  table  23.  With  all  quantities 
measured  in  angstrom  units,  S  =  co£( aizi/Rf )  ,  where  C0  *  116,140  and 
the  units  o^  S  are  cm-1 /(A)3.  The  polarizability  of  flourine  is  assumed 
to  be  1.0  A  .  Experimental  values  of  a2  were  computed  from  the  results 
of  table  22  and  the  appropriate  value  of  S.  If  it  is  assumed  that  the 
values  of  are  dependent  on  the  ion  only,  these  a 2  should  be  rela¬ 
tively  host  independent.  With  this  assumption,  the  average  value  for  o2 
given  in  table  24  can  be  used  to  calculate  the  energy  shifts  of  any 
host,  if  the  x-ray  positions  and  the  polarizability  of  the  ligands  are 
known. 


TABLE  22.  EXPERIMENTAL  ENERGY  (IN  UNITS  OF  103  cm*1)  OF  LOWEST 
ENERGY  LEVEL  OF  4fN"15d  TOR  FREE  ION,  AND  Af(J 


Ion 

Free  lona 

LaF3b 

YF3b 

luf/ 

LiYF4b 

Nd 

70.1 

60.0 

58.4 

57.1 

54.6 

Er 

75.4 

64.7 

64,3 

63.8 

62.5 

Tm 

74  „  14 

64.4 

63.4 

63.0 

61.5 

^Sugur  and  Reader,  1973,  Chem.  Phys.  59,  2083 . 
bYang  and  Deluca,  197 6,  Appl  Phys.  Lett.  29,  <99. 
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TABLE  23.  LIGAND  DISTANCE,  R  (A),  AND  MULTIPLICITY,  Z,  FOR 
COMPOUNDS  LISTED  IN  TABLE  22a 


Compound 


Ligand  distances  (multiplicity) 


LiYF4 

2.245 

(x4) 

2.293 

(  x4) 

LaFj 

2.246 

(x2) 

2.416 

(x2> 

2.443 

(xl  ) 

2.49 

(x2) 

2.64 

(x2) 

YF3 

2.253 

(xl) 

2.253 

(xl  ) 

2.266 

<x2> 

2.299 

(x2) 

2.323 

(  x2) 

2.595 

(xl  ) 

LuF-j 

2.191 

(xl) 

2.243 

<x2) 

2.257 

(x2> 

2.257 

(xl  ) 

2.277 

(x2) 

2.618 

(xl  ) 

aX-ray  data  used  in  these  calculations  are  from  C.  Keller  and  H.  Smutz, 
1965,  J,  Inorg.  N ucl ,  Chem.  27_,  900,  for  LiYF^;  A.  Zalkin,  D,  H.  Templeton, 
and  T.  E.  Hopkins,  1966,  Inorg .  Chem.  5_,  1466,  for  LaF^/  and  A.  Zelkin  and  D. 
H.  Templeton,  1953,  .7.  Am.  Chem.  Soc.  7j_,  2453,  for  YFj  and  LuFj. 


TABLE  24.  SUM  S  (cm'VA2)  OF  EQUATION  (313)  FOR  VALUES  IN 
TABLE  23,  Aj?d  -  Afd  FROM  TABLE  22,  AND  o 2  COMPUTED  FROM  EQUATION  (311) 


Compound  S 

Nda 

Er* 

Tm* 

Afd-Afd  °2 

A?d‘‘A£d  °2 

^d-Afd 

°2 

LiYF4 

6815 

15500 

2.274 

1  2900 

1  .893 

12800 

1.878 

LaF-j 

5184 

10100 

1 .948 

10700 

2.064 

9900 

1 .910 

YF3 

6898 

1 1  700 

1 .696 

11100 

1  .609 

10900 

1 .580 

LuF3 

7516 

13000 

1 .730 

11600 

1.543 

1 1  300 

1.503 

av  a2 

- 

- 

1.912 

- 

1  .777 

- 

1 .718 

aYang  and  DeLuca,  1976,  Appl.  Phys.  Lett.  29,  499 . 
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The  results  obtained  for  the  three  ions  studied  here  agree 
qualitatively  with  what  would  be  expected  from  the  usual  lanthanide 
contraction  (P.eisfeld  and  Jorgensen,  1977),  that  is,  (^(Nd^)  >  (^(Er^*) 
>  02(Tm3+)«  A  more  extensive  comparison  of  equation  (310)  with  the 
experimental  data  in  oxide  and  chloride  host  materials  is  necessary 
before  any  quantitative  claims  can  be  made.  Nevertheless,  the  results 
indicate  that  equation  (310)  can  be  a  useful  rule  of  thumb  indicating 
the  host  dependence  of  the  energy  spacings  between  the  various  configu¬ 
rations  of  rare-earth  ions. 


The  results  here  can  be  used  to  predict  a  host-dependent  spin- 
orbit  parameter  by  using  equation  (309)  in  the  usual  spin-orbit  Hamil¬ 
tonian  for  a  central  potential: 


,2  2  Z<°h 
h  S  V  1  i 

2  2  f  6 

m  c  l  R, 
i 


(312) 


Similarly,  the  result  given  can  be  used  to  give  host-dependent  Slater 
parameters  as 


AF 


(k) 


„  2 
aiZie 
R2k+4 
1 


(313) 


where  the  various  quantities  are  the  same  as  in  equation  (311).  The 
result  given  in  equation  (313)  has  been  noted  (Copeland  et  al,  1978)  for 
a  general  two-electron  interaction  and  is  intimately  related  to  a  pre¬ 
vious  result  where  the  ligands  are  replaced  by  an  isotropic  solid.  In 
view  of  the  difficulty  of  choosing  the  parameters  used  in  the  result 
given  ir  Morrison  et  al  (1967),  it  would  appear  that  equation  (313)  is  a 
preferable  form. 
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The  following  equations  include  the  most  commonly  used  relation¬ 
ships  needed  for  arriving  at  the  results  given  in  the  main  body  of  the 
report. 

A-1.  CLEBSCH-GORDAN  COEFFICIENTS 


A-1 . 1  Relation  to  3-j  Symbol 


<a(a)b(0)|c(Y)>  -  /5c+T  (-1)“a+b_Y(®  g  (A-1) 


A-1. 2  Symmetry 


a+b’c<a(-a)b(-0)|c(-Y)> 

(A-2) 

a+b“C<b(g)a(0) | c ( y ) > 

(A-3) 

<a(a)c(-Y)|b(-0)> 

(A-4) 

-  (-1}  <*wct-y 

(In  eq  (A-1),  (A-2),  (A-3),  and  (A-4),  a  +  0  «  y. ) 
A-1. 3  Orthogonality 


6__,  “I  <a(a)b(y-a)  |c(y)xa(a)b(t-a)  |c' (y)>  (A-5) 

cc  a 

4  ,i  ,  *1  <a(a)b(y-a) |c(Y)><a(a’ JbtY'-a* ) |c(y* )>  (A-6) 

aa  yy  c  ' 


A-1 . 4  Special  Values 

<a( 0)b( 0)  |c(0) >  -  (-1)8/2  (-gf-f4)l/2  “7 — w^w — T 

8  t(3i)t(s2)t(s3) 


(A-7) 


where 


»  -a+b+c,  s,  =  a-b+c,  s,  *  a+b-c,  s  *  a+b+c 


T(s)  a  ~r~ 

/ST 

<a(a)b(0) ja+b(o+0)> 


Y  2a  V 

Va+o/Vb 


/  2a+2b  \ 
(a+b+a+0) 


( A-8) 
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<i ( 0)k+2( 0) |£(0)> 


■awMaai 1/2 


<H (0)k ( 0)  Jt(0)> 


(A-9) 


WIGNER-ECHART  THEOREM 


A-2.1  Tensor  Tk<?  (in  L  Space)  Only 


<L'M'  |Tkg|LM>  =  <L(M)k(q)  |L»  (M' )><I.'llTkllL> 


(A-1 0) 


A-2.2  Mixed  Tensor  T, 


k  in  S  space  projection  X 
k  in  L  space  projection  q 

<I*'MLS'MslTXqlLMLSMS> 

=  <l( MjJ k (q)  I L'( M£)  ><s(Mg)x(X  )|s'(m|)> 
x  <L'S'IITK:kllLS> 


(A-1 1 ) 


A-2. 3  Matrix  Elements  of  Crystal  Field 


CJ'M'L’S'  |  H  |.TMLS> 


=  1  B*  <J(M)k(q)  |  J' (M' )><J'L'S#  U(K)II  JLS>/ 2£+1  <A#C  IIJI>  (A-1  2) 

kq  Kq  K 


CJ'L'SU  U(k)ll  JLS>  -  /  2J+1W(kLJ'S;L'J)  (L'Sl  LS) 


(A-1 3) 


(The  matrix  elements  (L'SU  LS)  are  tabulated  in  Neilson 
and  Koster,  1964.) 

A- 2 . 4  Matrix  Elements  of  Spin-Orbit  Interaction 


<J'M'L’S' |  K2|  JMLS>  =  [JUJl  +  1  )  (2£  +  1  )  l1/2 


x  W(1SL'J;S'L) 
x  (L'S'IV11ULS) 


(A-1 4) 
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(The  matrix  elements  (L'S'II  V* 1 8  LS)  are  tabulated  in 
and  Roster,  1964.) 

A-3.  COMMUTATION  RELATIONS 

A-3.1  Angular  Momentum  Operators 

[*M(i),Av(j)]  = /' 2  <1  (v)1  (u  )  1 1  (U+v)>£|J+v(i)6i;. 

[sydhSytj)]  =  /T  <1  (v  )1  (p  )  1 1  (U+v)>sy+v(i)6ij 

1  =  l  t( i),  UI  S(i>,  $  =  t+S 

i  i 

=  [yj>^J  -  [yjj.yj)] 

[  SP  (  j  }  »  Jy]  -  [8y(j)*Sv]  "  [S|j<j),Sv<j)] 

A-3. 2  Spherical  Tensors 

=  C^rjJ  *  yj 2j^+i  Ykq(ri) 

[  ,  C^q  ( i  )  ]  »  /k(k+1  )  <k(q>1(|l)|k(q-Hi)>C^  (i) 

[Vy,f(r)Ckq]  <1 (u )k(q) |  k 1 (q+y )><1 (0)k(0) |k' (0) >Ck , ^ q+^ 

*  [ir-  +  ^  f (r)[k(k+1  )(2k+1  )]1/2  W(klkMjkl) 

[V^r^q]  =  -[k(2k+1  )]1/2  <k(q)1(U)|k-l(q4|i>>rX”1Ck_1 

[V),,r-|'-'c)lq] 

a  -[  (k+1 )  (2k+1 )] 1  /2  <k(q)1(y  )|k+1(q-Hi)>r_k‘2Ck+1fq+iJ 


Neilson 

(A-15) 

( A-1 6) 

( A-1  7 ) 

( A-1 8 ) 

( A-1 9 ) 

( A-20 ) 

( A-21 ) 

( A-22 ) 

( A-23 ) 


In  equations  (A-21)  through  (A-23),  C^g  *  Ckq(r). 
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A-4.  RACAH  COEFFICIENTS  AND  9-j  SYMBOLS 
A-4.1  Working  Definition 

<a(a)b( 6) | e( a+6 ) ><«( a+0)d( <5 ) |c(a+8+6)> 


l  /T3e+T H 27+TJ  W(abcd,*ef  )<b(  0 }d < 5 )  |  f  ( 0+6 ) ><a(a)f  ( 0+5)  |c(  a+8+6 ) > 


(A-24) 


W( abed; Of)  ■ 


/  (  2b+1  )  (  2d+1 ) 

A-4. 2  Orthogonality 

£  (2e+1 )(2f+1 )W(abcdjef )W(abcd;eg)  *  6^ 
e  " 

A-4. 3  Relation  to  6-j  Symbol 


W(abcd;ef )  »  (-1 ) 


A-4. 4  9-j  or  X  Symbol 


a+b+c+d  /a  b  e 
Id  c  V 


X(abc,def ,ghi )  =»  £  (2k+1 )W(aidh>kg)W(bfhd>ke)W(aibf ;kc) 
k 

6  -6  .(-1 )C+g“a_eW(abde jcg) 

X(abc,def ,ghO)  = - ° - - 

/(2c+1 )<2g+1 ) 

A-5.  SPHERICAL  BASIS  VECTORS 
A-5.1  Unit  Vectors 

e±1  -  (©x  ±  iSy)//I 

a 

e  »  e 
o  z 

eS  -  (-1)Pe-iJ 

A  A 

6**e  »  6 

U  v  uyv 


eiJxev  =-  -1/7  <1  (v)1  (y)  1 1  (y+v)>e 


y+v 


(A-25) 


(A-26) 


(A-27) 


(A-28) 


( A-29 ) 


( A-30) 
( A-31 ) 
(A-32) 
(A-33) 
(A-34) 
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A-5.2  Examples 

*  - 1  Hk 


( A-35 ) 


t  -  y 


(A-36) 


*•*  *  I  *2%  “  £ 


A  ^  A  A 

r  -  J  ^5ciu<r> 


A-5 . 3  Generalizations 

Lot  A  ,  B.  etc,  be  spherical  tensors  and  define 
«s  bp 

A  "B  -  y  (-1 )°  A  B  «  B  »A  . 

*•  act  a,— a  ^-a  ~a 


(A-37) 


(A-38) 


(A  and  B  commute. ) 


{^a  SbU  *  l  <a(a)b(Y-«)!c(Y)>  VW* 


VRa  -  (-l)«  /5a?T 

teaSblc’Sc  *  <’1>a+b"CMa}c*Sc 

*  UaSclb-Bb 

“  < -1  )a  ViSr  {£cfib}*&a  * 


A-6.  EXPANSIONS 


A-6.1  Addition  Theorem 


Pk(r*R)  -  J  Ckq(^)Ctq(ft) 

(P^  is  a  Legendre  polynomial.) 

Caa(r)Ct>g(r)  ■  |  <a(0)b(0)  |k(O)xa(a)b(0)  Ikta+BJX^ >cj+g(r) 


(A-39) 


( A-40) 
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Wr1,r2) 


-  I  <a(«)b(tHO|K<e»Ca#0(r1)<^r,H.(r2) 

-  [^fr1)^b(r2)]KQ 


(A-41) 


A-6. 2  Sinqle-Centei  Expansion 


i_  v  *  ckg(r)C£g(R> 


|*4 


Tc+r 


(A-42 ) 


k  R 


kTTCk(r)‘Ck(R) 


(A-43) 


A-6. 3  Two-Center  Expansions 


J _  .  V  f  2a+2b-)1/2  r c  (*)C  ,".1  #c  fR)  _x&yb 

-x-t|  5bl  2a  a'x '  Dly  Ja+b  a+blR)  Ra+b+1 


(A-44) 


|*44|  "  lq  |*-i|k+1  Ckq(y)Ci?q(R’X) 


(A-45) 


Ckq(R-X) 

|*-x|k+1 


"L  (22a2k)l/2  <a(«^(q)U+K(a+q)>C*a(x)Ca+kfa+q(R) 


(A-46) 


aa 


A-6. 4  Miscellaneous  Expansions 

la  ( 2a+1  )ia(Xr)ka(XR)Cao(r )C*a(R)  (A-47) 


(in(z)  =  J  "ft  f  2z  In+i  /2(z *n(z)  m  ^  /^z  ^n+l /2^z  ^ ;  ^  and  K  are 
Bessel  functions.)  '  * 


iX  |*~r I 

S--'j|3|  |  -  j.  (2a+1)ja(Xr)[ya(XR)  +  ija(XR)]caa(^ )C*Q(R)  (A-48) 


(j.  and  ym  are  spherical  Bessel  functions.) 

el  3 
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APPENDIX  B 


8-1 .  INTRODUCTION 

The  dominant  quantity  in  electrostatics  is  the  electric  poten¬ 
tial.  Once  having  found  the  potential,  we  can  obtain  all  other  essen¬ 
tials  by  suitable  operation  on  the  potential.  In  particular,  if  the 
potential,  ${r),  is  given,  then  the  electric  field  is  obtained  by  taking 
the  gradient  of  $(r),  that  is, 

I  =  -V$(r)  .  ( B-1 ) 

In  equation  (B-1),  it  is  assumed  that  <j>  is  a  function  of  the  general 
field  coordinates  r>  however,  since  Q  may  also  involve  other  coordinates 
in  a  complicated  manner,  care  must  be  taken  to  separate  the  field  point 
from  other  coordinates  so  that  the  gradient  in  equation  (B-1)  i9  with 
respect  to  r.  We  shall  always  here  interpret  the  gradient^  operator  in 
equation  (B-1)  to  be  taken  with  respect  to  the  field  point  r. 


B-2.  EXAMPLES 

A  few  examples  will  suffice  to  clearly  illustrate  the  above 
points. 


B-2.1  Point  Charge  at  Origin 

If  there  is  point  charge  q  at  the^  origin  of  our  coordinate 
systom,  then  the  potential  at  the  field  point  r  is  given  by 


♦ 


l?l 


r 


(xz  +  y2  +  z2)l/2  .  (B-2) 

The  electric  field  for  this  simple  system  is  given  by  -grad  4),  and  is 

•* 

E  -  -J^-r  ,  (B-3) 

l?l3 

since 


1 


♦ 

r 


(B-4) 


1*1  r3 

The  relation  given  in  equation  (B-4)  can  be  worked  out  in  detail  by 
taking  each  vector  component  of  V  as 


l?l 


*«  fx  -jjf  +  ey  fa  -[jJ  +  ez  fa  -j7 


where  jr|  is  given  by  aquation  (B-2). 
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The  potential  of  a  dipole  at  the  origin  can  be  obtained  by 
a  slight  extension  of  equation  (B-2).  If  we  assume  the  dipole  is  made 
up  of  two  oppositely  charged  particles,  q  and  -q,  located  symmetrically 
about  the  origin  at  1/2  and  -t/2,  respectively,  the  resulting  potential 
using  equation  (B-2)  is 


I  r-e/2 i 


I  r+e/2 | 


( B-6 ) 


By  noting  that 


j?+£/2|  *  [  (r+e/2)»  [i+t/2)  ]  =  (r+e»?)1//2 


(f2+?.?U/2  « 

2r' 


and  substituting  these  results  into  equation  (B-6),  we  get 


( B-7 ) 


It  is  customary  to  let  e  0  and  q  +j*>  in  such  a  way  that  qe  remains 
finite  and  is  the  point  dipole  moment  p.  Then  equation  (B-7)  gives 


♦  ,  ( B-8) 

r3 

which  can  be  and  is  taken  as  the  potential  at  the  field  point  r  due  to  a 
dipole  of  moment  p  at  the  origin.  The  potential  given  ir,  equation  (B-8) 
can  be  written  in  a  more  compai  t  and  frequently  more  useful  form  using 
equation  (B-4)  or 

#d  -  -p»7  -J—  .  (B-9) 


B-2 . 2  Point  Charge  ai  R 

We  consider  a  point  charge  q  now  located  at  a  point  from 
the  origin;  we  wish  to  evaluate  the  potential  at  the  field  point  r. 
Then 


I*-?  I 


( B-1 0 ) 
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B-3.  ATOMS  AND  MOLECULES 


We  now  wish  to  find  the  energy  of  interaction  of  two  charges  in 
the  presence  of  a  remot^  polarizable^ion.  The  coordinate  system  is  such 
that  the  charges  q.,  at  r^  and  q2  at  r.  have  a  common  origin.  The  polar¬ 
izable  ion  of  charge  Q  is  located  at  £  such  that  R  >>  r1  or  r,*  The 
electrostatic  potential  at  r  due  to  the  charge  ,  given  by  equation  (3- 
2),  is 

q 

kj-rl 

We  use  two  arguments  in  +the  potential  to  signify  the  source  point, 
r^ ,  and  the  field  point,  r.  If  there  were  no  charges  or  sources  other 
than  q1  and  q2,  the  energy  of  interaction  would  be 


»(M2) 


( B-1 7 ) 


a  familiar  result.  Thus  we  can  write  the  energy  of  a  point  charge,  q^( 
in  a  field  whose  sources  are  independent  of  q^  as 


°(?i)  -  v(?)  -  ♦  • 

r=r . 


(B-1 8) 


With  these  results  the  energy  of  the  system  due  to  the  point  charges 
only  is 


°  IV?2I 


( B— 1 9) 


The  energy  given  in  equation  (B-1 9)  does  not  include  the  energy  due  to 
the  polarization  of  the  ion  at  $.  To  evaluate  the  polarization  energy 
we  must  first  find  the  dipole  moment  induced  in  the  ion.  From  equation 
(B-12)  (with  &  *  r1 )  we  have  for  the  potential  at  r  due  to  q^ : 


mTTT7  ' 

r~ril 


(B-20) 


and  from  equation  (B-1 3): 
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E(r) 


+  4 

r-j* 


1  I  ♦  +  I  3 

lr-r1 1 


The  dipole  moment  at  &  in  given  by 


4- 

P 


aE(r ) 


=  aq 


1 


3 


(B-21  ) 


( B-22 ) 


where  3  is  the  polarizability  of  the  ion  at  R. 

Hie  potential  at  r  due  to  this  dipole  at  £  is  given  by  equation 
( B- I  4 ) : 

4>d(r)  «  -p  •  ( B-23 ) 

I*-?!3 


and  the  energy  of  interaction  with  the  charge  q2  is 

ud(V2)  -,,♦<*>! 


4- 

I  r=r_ 


(B-24) 


from  equation  (B-17).  Substituting  equations  (B-23)  and  (B-22)  into 
equation  (B-24),  we  have 

(?'  2) '  ‘Wl’q2  ' 


(B-25) 


which  is  symmetric  in  the  two  charges  and  positions.  Thus,  we  could 
have  obtained  the  identical  renult  if  we  had  used  charge  q 2  and  its 
field  to  find  the  induced  dipole  and  evaluated  the  potential  due  to  this 
dipole  at  r^ .  It  is  possible  to  extend  these  results  to  more  compli¬ 
cated  systems,  but  if  this  is  seriously  intended  the  more  elegant  tech¬ 
niques  given  by  Judd1  should  be  used.  The  result  given  in  equation  (B~ 
25)  was  used  to  develop  a  two-electron  crystal-field  interaction2 
following  a  suggestion  of  Judd. 


1B.  R.  Judd,  Nath.  Proc.  Camb.  Phil.  Soc .,  80_  535  (1976). 

2C.  A.  Morrison,  J.  Chem.  Phys.  72,  1001  (1930). 
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B-4.  POSSIBLE  APPLICATION  TC  CATALYSIS 


Perhaps  a  more  important  application  of  equation  (B-25)  is  the 
application  of  this  result  to  the  theory  of  homogeneous  catalysis.  In 
homogeneous  catalysis  an  inert  gas  (Q  *  0  in  eq  (B-19))  is  introduced 
into  a  mixture  of  gases  represented  by  q1  and  qj.  We  assume  here  that 
q1  and  q£  represent  the  charges  on  atoms  1  and  2  and  the  inert  gas  is  to 
speed  up  a  desired  compound  formed  by  some  combination  of  q1  and  q2* 
The  interaction  energy  given  by  equations  (B-25)  and  < B— 17)  is 


q,q2 


M1q2 


{*-?,) -(H) 

|S-?,|3|S-f2|3 


( B-26) 


We  see  that  if  q1  and  q2  are  of  the  same  sign,  the  interaction  with  the 
catalyst  at  R  is  such  as  to  reduce  the  repulsion.  To  apply  equation  (B- 
26)  to  the  present  case  we  can  take  r *  with  no  loss  in  generality, 
and  equation  (B-26)  becomes 


,„i  i,  qlq2  _  <$-r/2).($+r/2) 

U  ( r ,  R )  =  I  r  -  aq  q  — — - - 7  • 

ir'  1  2  I i£-r /2  I  3  I  $+r /2  !  3 


(B-27) 


with  r  =  r 


>  _> 
C  *  2V 


|S-?/2|  3|5«/2| 
where  r  is  the  distance  between  atoms  1  and  2. 


Thus  if  the  two  atoms  have  the  same  sign,  the  repulsive  forces  are 
reduced;  this  reduction  may  be  sufficient  to  allow  the  two  atoms  to  get 
close  enough  to  permit  electron  exchange,  which  may  be  necessary  to 
speed  up  a  desired  compound  formed  by  some  combination  of  q^  and 


B-5.  ELECTROSTATIC  POTENTIAL  DUE  TO  MULTIPOLAR  DISTRIBUTION 

In  much  of  our  work  we  shall  need  the  expression  for  the  electro¬ 
static  potential  due  to  a  distribution  of  multipolar  moments.  A 
straightforward  derivation  of  the  desired  results  can  be  obtained  by  an 
application  of  the  simple  laws  of  electrostatics  given  above.  If  we 
have  a  charge  distribution,  p(x),  at  the  origin  of  a  coordinate  system, 
then  the  electric  potential  at  5  from  the  origin  is 

p(x)dx 

d<fr(r)  = — — —  ,  (B-28) 

I  ♦  “*■  \ 
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where  dT  Is  the  volume  element  at  x  from  the  origin.  The  denominator 
in  equation  (B-28)  can  be  expanded  in  Legendre  polynomials  as 


— - —  =■  y  — —  p  r  cos  9  ) 

,  L  .  n  n1-  xrJ 

|  r-x  j  n*0  r 


(B-29) 


where  0^  is  the  angle  between  x  and  r  and  we  have  assumed  |r|  >>  |£| 

By  Legendre's  addition  theorem, 


P  f  cos  0  } 

n^  xrJ 


n 

!  C  (r)C*  (r)  , 

nm  nm 

m=-n 


( B-30) 


A  A  ^ 

where  x  and  r  are  unit  vectors  along  x  and  r,  respectively,  and  x*r  ■» 
cos(0rx).  If  we  substitute  equation  (B-30)  into  equation  (B-29),  we 
have 


„  n  a  a 

— ! =  V  c  (x)C*  (r)  . 

i+  ♦  i  L  n+1  nm  nm 


(B-31 ) 


I  r-x  nm  r 


In  equation  (B-31 )  we  have  dropped  the  limits  on  the  sum,  a  technique  we 
shall  use  throughout  unless  the  sums  are  restricted,  in  which  case  the 
limits  shall  be  explicitly  given.  If  the  result  given  in  equation  (B- 
31)  is  substituted  into  equation  (B-28)  we  have 


•Mr)  »  /  dT  p  (x)  l  C(x)C*  (r) 

x  n+i  nm  nm 

n,m  r 


( B-32 ) 


Integrating  the  results  in  equation  (B-32)  produces 


♦  <t>  «  I  CJm(r,  , 

n.m  r 


( B-33 ) 


where  we  have  defined  the  multipolar  moment  as 


0  «  /  p(x)xnC  (x)  dT 

Tim  nm  i 


(B-34) 


The  definition  of  the  nth  multipolar  moment  given  in  equation  (B-34) 
agrees  with  the  customary  definition  for  n  *  1  (dipole  moment)  but 
differs  from  a  number  of  others  for  .1  =  2.  Frequently  Q2m  is  defined  by 


82m  "  I  •>•*>[  2* 2<V>»J  dTx  ’ 
we  shall  use  only  the  definition  given  in  equation  (B-34). 


( B-35 ) 
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Since  the  electric  potential  is  real,  equation  (B-33)  can  be 
written 


<j>  (r ) 


I 

nm 


Q*  C  (r) 
nm  nm 

n+i 


(B-36) 


We  shall  use  the  form  given  by  equation  (B-33)  or  equation  (B-36)  de¬ 
pending  on  which  is  the  more  convenient  in  a  particular  problem. 

If  the  charge  distribution  p(x)  is  located  at  r  from  the  origin, 
the  potential  is  given  by 


<J>(r)  *  / 


P(x)  dt 


|r+xj 


<K  r ) 


I 

nm 


(-1  )  Q*  C  (r) 
_ nm  r.m 

n+1 


(B-37) 


(B-38) 


a  result  which  could  have  been  obtained  from  equation  (B-36)  by  letting 
r  ♦  -r  and  noting  that  Cnm(-r)  =  (-1)nCnm(r).  It  should  be  noted  that 
the  replacement  r  >  -r  in  equation  (B-28)  converts  the  denominator  there 
into  the  denominator  in  equation  (B-37). 

We  then  have  from  equation  (B-36)  the  electric  potential  at  r  due 
to  a  charge  distribution  at  the  origin,  and  from  equation  (B-38),  we 
havj  the  electric  potential  at  the  origin  due  to  a  charge  distribution 
at  r. 


A  second  and  independent  definition  of  Qnm  is  given  by 

0  -  «  E  ,  ( B-39 ) 

him  n  nm 

where  a  is  the  nth-pole  polarizability  and  Enm  is  the  nth  multipole- 
inducing  field.  For  n  =  1  (dipole)  we  have 

-Vi.  ■  <B-40) 

which  is  the  spherical  representation  of  the  relative 

p  *  all  ,  (B-41 ) 

where  a  »  a  is  the  ordinary  polarizability  and  E  the  ordinary  electric 
field.  If  we  write  the  vector  relations  in  equation  (B-41)  in  spherical 
form,  we  have 
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Similarly,  if  we  have  the  charge  qQ  at  the  origin  and  wish  the  multi¬ 
polar  inducing  field  at  5,  the  multipolar  inducing  field  is 


nm 


=  -q 


,n  +  1 


( B-49) 


If  we  have  an  ion  at  the  origin  whose  nth-pole  multipolarizability 
is  a^,  the  induced  moment  Qnm  is  given  by 


Q  ■-  a  E  , 
xnm  n  nm 

and  from  equation  (B-48)  we  have 

C  (R) 

_  nm 

®nm  ~  an^0  Rn+1 


( B-50) 


( B-51 ) 


The  potential  at  the  point  r  of  this  induced  dipole  at  the  origin  is 
given  by  equation  (B-36)  to  obtain 


-  C*  (R)  C  ( r ) 

4>(r)  =  nm  "  ~3+T- 


( B-52) 


where,  to  repeat,  5  is  the  position  of  the  charge  qQ  and  r  is  the  field 
point. 


The  following  +  problem  illustrates  all  of  the  above.  Assume  a 
point  charge  qQ  at  r  and  a  dielectric  sphere  of  radius  a  at  the  origin 
with  dielectric  constant  e.  The  potential  in  the  two  regions  is 


and 


4>,  (r)  =  l  a*  rnC  (r)  ,  r  <  a  , 

l  u  nm  nm 


*o(r) 


b*  C  (r) 
nm  nm 

n+1 


r  <  a  . 


By  equating  the  potential  at  r  =  a  and  the  derivatives 


e 


*3r~ 


r=a 


3*o 

8r 


r=a 
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(normal  components  of  D  are  continuous),  a_„  and  can  be  determined. 

nm  nm 

From  equation  (B-52)  we  have 

h  .  .  anqOCnm(g) 

“nm  j^n+1  ' 

if  we  take 


an  =  [n(e  -  1)a2n+1]/(en  +  n  1 ) 

Notice  that  this  an  agrees  with  tha  an  given  by  Judd1  (Judd's  K  is 
our  e ) . 


If  we  assume  that  and  a  are  known,  then  we  can  eliminate  e 

from  a  to  obtain 

n  _  2n+i 

Snc^a 

a  *  - > 

(n-1)«1  +■  (2n+1)a 


which  could  be  used  to  estimate  the  a  for  n  >  1  for  an  ion  whose  a1  is 
known  and  a  might  be  taken  as  the  ionic  radius.  Comparing  equation  (B- 
52)  with  the  expression  for  <}>o  we  have 


b 

nm 


-a  q„C  (R)/R 
n  0  nm  ' 


n+1 


1B.  R.  Judd,  Math.  Proc .  Camb.  Phil.  Soc.,  80  535  (1976). 
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